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ABSTRACT 

We use the overdensity field reconstructed in the volume of the COSMOS area to study the nonlinear 
biasing of the zCOSMOS galaxies. The galaxy overdensity field is reconstructed using the current 
sample of ~8500 accurate zCOSMOS redshifts at Iab < 22.5 out to z^l on scales R from 8 to 12 
h~^Mpc. By comparing the probability distribution function (PDF) of galaxy density contrast Sg to 
the lognormal approximation of the PDF of the mass density contrast i5, we obtain the mean biasing 
function b(6, z, R) between the galaxy and matter overdensity field and its second moments b and b up 
to z ~ 1. Over the redshift interval 0.4 < z < 1 the conditional mean function {Sg\S) = b{S, z, R)S is 
of the following characteristic shape. The function vanishes in the most underdense regions and then 
sharply rises in a nonlinear way towards the mean densities. {Sg\6) is almost a linear tracer of the 
matter in the overdense regions, up to the most overdense regions in which it is nonlinear again and the 
local effective slope of {Sg\S) vs. S is smaller than unity. The {Sg\S) function is evolving only slightly 

over the redshift interval 0.4 < z < 1. The linear biasing parameter increases from b = 1.24 ± 0.11 
a.t z = 0.4 to 6 = 1.64 ± 0.15 at z = 1 for the Mb < —20 — z sample of galaxies, b does not show 
any dependence on the smoothing scale from 8 to 12 h~^Mpc, but increases with luminosity. The 
measured nonlinearity parameter b/b is of the order of a few percent (but it can be consistent with 0) 
and it does not change with redshift, the smoothing scale or the luminosity. By matching the linear 
bias of galaxies to the halo bias, we infer that the Mb < —20 — z galaxies reside in dark matter haloes 
with a characteristic mass of about 3 — 6 x 10^^ Mq, depending on the halo bias fit. 
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1. INTRODUCTION 

The large scale structure in the Universe is believed 
to have formed via gravitational instability of small, 
primordial density fluctuations. Virialised dark matter 
haloes are produced by the collapse of some overdense 
regions followed by hierarchical merging. Galaxies are 
formed within the dark matter haloes through multi- 
plex processes including gas cooling, star formation and 
feedback which are difhcult to model accurately (e.g. 
I White k Reeslll978f) . It is thus expected that the re- 
lation between galaxies and the underlying matter dis- 
tribution will be also complex. In particular, the effi- 
ciency of galaxy formation and the rate of galaxy evo- 
lution might vary from place to place depending on the 
matter density field. Therefore the actual galaxy dis- 
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tributions is not expected to be a rightful tracer of the 
underlying mass. This phenomenon is often referred to 
as "galaxy biasing" . 

In order to extract cosmological information from 
galaxy surveys it is important to model and parame- 
terise galaxy biasing. This often requires using a sta- 
tistical approach. Assuming that galaxies preferentially 
form within the pe aks of the primordial density distribu- 
tion, [K^se^ (|1984[ ) showed that the two-point correlation 
function of the galaxy distribution, ^gg, should be ampli- 
fied with respect to the mass autocorrelation function, 
^mm, according to the relation: 



and 



(r) 



(1) 



where the "biasing parameter" is independent of the 
spatial separation r. A similar relation is obtained by 
relating the density contrast of the matter 6 and of the 
galaxies Sg at some position r through the deterministic 
and linear relation 



b6{r) 



(2) 



This is the simplest model for galaxy biasing and is still 
commonly used. While equation [T] follows from equa- 
tion[2l the opposite does not hold. An obvious deficiency 
in the definition of dg above is that it will break down in 
the most underdense regions (5 ^ if 6 > 1, as values of 
(5g < — 1 are not possible. This implies that the galaxy 
bias b must be a nonlinear function of 5 and, in gen- 
eral, it can also vary with redshift z, galaxy type and the 
smoothing scale R used to define the density contrast: 

b = biS,z,R). (3) 

iFrv fc: Gaztanagal ()1993l ) proposed to parameterisc this 
function in terms of coefficients of the Taylor expansion 

Sg = bf, + bi6+^6^ + ... , (4) 

which are not fully independent as the conditions (Sg) = 
and Sg{d = —1) = —1 must hold. 

Galaxy biasing is also expected to have a stochastic 
element: for any given value of 5 there will be a whole 
distribution of values for Sg. The stochasticity originates 
from a number of different sources. First, the dynamics 
of large scale flows depends on extra variables beyond 
the value of the local density contrast 6 (e.g. on the 
tidal tensor) and mak es the bias relation n onlinear, non- 
local and stochastic (jCatelan et al.|[r998[ ). Second, the 
efficiency of galaxy formation depends on details of the 
gas physics. Third, galaxies are discrete objects and any 
attem pt to recons truct Sg will be effected by shot noise. 

De kel fc Lahavl 11999) have proposed a formalism 
which separately accounts for the nonlinearity and 
stochasticity of the biasing process. Galaxy biasing is 
described in terms of the conditional probability func- 
tion P{5g\S) and its moments. A key quantity here is 
the mean biasing function b{d) defined by the conditional 
mean: 



b{S)S=(Sg\S) = / dSgP{Sg\S)S.^ 



(5) 



The mean biasing function b{6) and its nonlinearity can 
be characterised by its second non-trivial moments: 



(7) 



with cr^ f}^f2 variance of the mass density contrast distri- 
bution. The parameter b measures the slope of the linear 
regression of 6g against S. In the case of linear biasing 

(see equation [5]) , both b and b reduce to the constant 

bias. The ratio b/b is thus a measure of the nonlinearity 
in the biasing relation. Moreover, the local variance of 
5g at fixed S, crg(i5) can be used to quantify the degree of 
stochasticity of the biasing relation. 

Based on the Press-Sche c hter formalism and its ex- 
tensions (|Bond et al.lll991h . iMo fc Whit^ (|1996l ) devel- 
oped an analytical model for the mean biasing relation 
of the dark matter haloes. This assumes that large 
scale motions follow the spherical collap se approxima 
tion. The general case is discussed by ICatelan et al 
(U998). Related work has been presented in iMo et al 
(|1997l) and Porciani et al. (1998, see also Scanna- 
pieco fc Barkana 2002) where two-point and higher- 
order statistics are considered. Follo wing the ana- 
lytical approach by iMo fc Whitd (|1996l ). a number of 
studies based on N-body sir nulations were carried out 
to st u dy the halo bias (e.g. Jin^ Il998t iPorciani et al.l 
T99I IShetii fc LemsonI 119991: ISheth fc TormenI 11999: 



Jhia 119991: IKravtsov fc KlvpinI [l 999: Shcth et aL 200i 
Seliak fc Warrenll2004l : lTinker et al...2005i : ,Pillepich et al 



{b{S)S^ 



(6) 



20081 ) ■ leading to a new set of fitting formulae for the 
mean biasing relation, and a better understanding of the 
origin of halo biasing. Independently of the exact halo 
definition, assumed cosmology, simulation box size and 
resolution, there is a consensus that in a cold-dark mat- 
ter scenario: i) at a given epoch, more massive haloes are 
more biased tracers of the underlying matter than lower 
mass haloes; ii) for halos of fixed mass, the amount of 
biasing increases with redshift. 

However, it is still a huge step from a successful de- 
scription of "halo biasing" to that of "galaxy biasing" 
as the latter requires incorporating a recipe for galaxy 
formation (and evolution) within the current cosmo- 
logical framework. Galaxy biasing has been studied 
through hydrodynamic simulations (e.g. iBlanton et al.l 
[19991 120001 : fGen fc OstrikeHl200fll : lYoshikawa et al EM) 
and semi-analytical modelling combined wit h N-body 
simulations (e.g. p<!a uffmann et al. 1997; B enson et al.l 
I2OOOI: ISomerville et alJfeOOlHSigad et al.ll2000i . hereafter 
SBD). Despite the difference in the treatment of the var- 
ious gas-related processes, all these studies reach the fol- 
lowing consistent conclusions: galaxy biasing is expected 
to be nonlinear, to depend on the properties of the con- 
sidered galaxy sample (and of their host dark matter 
haloes), and to be a function of cosmic time. 

There is now lot of observational evidence for galaxy bi- 
asing, and its dependence on galaxy type or redshift . The 
Dressler's morphology-density relation (Dressier 1980f) 
- the observational evidence that early type galaxies 
are more abundant in dense regions than spiral galax- 
ies - is a textbook example for this. Building on this, 
one can summarise decades of observations in the lo- 
cal universe with a statement that bulge-dominated, red 
galaxies with mainly old stellar populations preferen- 
tially live in dense regions, and they are more strongly 
clustered than the disk-dominated, blue, young galax- 
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ies (e.g. iNorberg et all |2001| . l2002t IZehavi et al]|2005D . 

Consistently, deriving the bias parameter from cluster- 
ing studies (see e.g. Equation[T]) suggests that early-type 
galaxies have a higher bias than l ate- type galaxies at all 
luminosities (Norb erg et al.l !20021. On the other hand, 
Hl-selected galaxie s have some of the low est bias values of 
all known objects (jBasilakos et al.ir2007f ). At higher red- 
shifts, galaxies are more biased tracers of matter, with 
linear bias paramet ers of b ^ 1.48 a t 0.7 < z < 1.3 for 
Mb < -2 galaxies (ICoil et al.l l2006l. see also lPoUo et"all 
120061 and iMeneux et al.ll2009D . up to the highest biased 
samples of extremely red objects (EROs), Lyman Break 
galaxies (LBGs) and Lyman-a emitters. For example, 
the bias parameter of EROs at z = 1.2 is & = 2.7 
(jMoustakas fc Somervillil2002f) . LBGs at z ^ 3.8 and 
4.9 is 5 ~ 2.5 and 6 ~ 4, respectively ()Lee et al.l 
120061 ) and the bi as parameter of L yman-a emitters at 
z ~ 4.5 IS 6-3.7(11 ovac et al]|2007f ). However, inferring 
the exact redshift evolution of the biasing process from 
the observational data is not straightforward, because 
galaxy surveys typically sample d ifferent population s of 
galaxies at different redshifts (e.g. iKovac et al.|[2007D . 

A simple way to model the biasing of galaxies which 
has received a lot of attention recently, is through the 
halo occupation distribution (HOD) formalism (e.g. see 
review by Cooray & Sheth 2002, and references therein). 
This method splits the bias problem into two steps: i) 
N-body simulations are used to characterise the spatial 
distribution and the clustering properties of virialised 
dark matter halos as a function of their mass (and/or 
some other properties) for a given cosmology; ii) the 
galaxy distribution is described in terms of the proba- 
bility distribution that a halo of mass M hosts N galax- 
ies of a specified type (the HOD). The first N moments 
of the HOD can be measured by fitting observed N- 
point statistics. Using the two-point correlation function, 
this approach has been widely employed to estimate the 
mass of the host halos of galaxies an d qua sars at low 
and high redshift (e.g. iMagliocchett i & Por cianil l2003t 
Pprciani et al. 2004; Abazaiian et al. 2005; Z ehavi et al.l 
2005; Phlcps et al. 2006; Zheng et al. 2007) and also 
to derive the mass-to-light ratio of virialised c osmic 
structures (e.g. lYang et^al.! 120051 : [Tinker et all 120051 : 
Ivan den Bosch et al.ll2007[ ). 

It is however interesting to go beyond the measurement 
of a (possibly sca le-dependent) bias paramete r from two- 
point statistics. iTegmark fc BromlevI (|1999) presented 
evidence that the present-day galaxy biasing is nonlinear 
and stochastic, employing the galaxy clustering in Las 
Campanas Redshift Survey. Using the 2 degree Field 
Galaxy Redshift Survey (2dFGRS) data, some contra- 
dictory results on the nonlinea r natur e of g alaxy bias 
have been obtained. While Verde et al.l (|2002[ ) found no 
significan t evidence for non l ineari ty from the bispectrum 
analysis, iGaztanaga et al.l ()2005f ) detect non-vanishing 
quadratic corrections in three -poin t correlation function . 
Moreover, iWild et al.l (|200l and iConwav et all (|2005l ) 
using the count-in-cells analysis exclude the determinis- 
tic linear bias model in both the flux limited and vol- 
ume limited (luminosity complete) samples of galaxies 
and find evidence for stochasticity. Measurements of 
the three-point correlation function and counts-in-cells 
statistics for galaxy samples from the Sloan Digital Sky 
Survey (SDSS) suggest that galaxy biasing is nonlinear 



and f airly complex ( Kavo et all 12004 iNishimichi et al.l 
12001 ISwanson et al.ll2008fl . 

The mean nonlinear biasing function of a sample of 
galaxies can be constrained by combining counts-in-cells 
measurements with models for the probability distribu- 
tion function (PDF) of mass density fluc t uatio ns (SBD, 
ISzapudi fc PanI l2004f) . iMarinoni et"all (|2005f) applied 
this techni que to the flrst-epoch VIMOS VLT deep sur- 
vey (VVDS lLe Fevre et al.ll2005D . over the redshift range 
0.4-1.5 on a characteristic scale R from 5 to 10 h~^Mpc. 
They conclude that galaxy bias increases with redshift 
and is nonlinear in all redshift bins probed. In addition, 
brighter galaxies are more strongly biased than less lu- 
minous ones, as well as redder galaxies are more biased 
than blue ones independently of redshift. However, the 
area covered by the VVDS used for this study was rather 
small (0.4 x 0.4 deg^) and the results are likely affected 
by cosmic variance. 

It is therefore very important to cross-check their ro- 
bustness against richer datasets. With this spirit, in this 
paper we perform a similar analysis on a larger sample of 
high redshift galaxies. We make use of the flr st — 10000 
spect ra of the zCOSMOS redshift survey (jLillv et al.l 
l2009f ) to measure galaxy densities within top-hat spheres 
of radius 8-12 h~ ^Mpc in the redshift range 0.4-1 (see 
IKovac et al.ll2009l ). Assuming concordance cosmology, 
we then use the PDF of the galaxy density contrast to 
estimate the mean bias function b — b{6, z, R) between 
zCOSMOS galaxies and matter overdensities. Within 
the limits of the zCOSMOS survey, we explore the shape 
and the luminosity, scale and redshift dependence of the 
conditional mean function {5g\5). A careful analysis of 
random and systematic errors based on mock galaxy cat- 
alogues is also presented. Finally, the characteristic mass 
of the halos hosting zCOSMOS galaxies is inferred from 
the value of the linear bias parameter. The structure of 
this paper is as follows. We present the data used for 
the analysis and describe the method to reconstruct the 
density field in Section [2j In Section [3] we provide the de- 
tails of the method which we adopt to derive the biasing 
function. We present the results of the tests on the mock 
catalogues in Section [H Finally, we present the biasing 
analysis on the real data set in Section [5l comparison 
to the literature results in Section [HI and conclusions in 
Section El 

The assumed cosmology is specified by erg = 0.8, 
17,„,o = 0.25, r2A,o = 0.75 and = 70 km s^^Mpc^i. 
We present the obtained results scaled to the units of 
dimensionless h parameter, Hq = lOO/i km s~^ Mpc~^, 
except of magnitudes, which are calculated with h = 0.7. 

2. THE ZCOSMOS SURVEY 

2.1. The 10k sample 

zCOSMOS (ILillv et al.ll2007l 1200^ is a redshift sur- 
vey undertaken in the fie ld of the mu l tiwave length 
2 deg2 COSMOS survey (iScoville et al.l HoOTa), the 
largest contiguous mosaic ever obtained using HST/ACS 
(|Koekemoer et al.ll2007l) . The zCOSMOS-bright sam- 
ple is a purely flux limited part of the survey, selected 
at Iab < 22.5. It covers the central ~1.7 deg^ of the 
COSMOS field. This corresponds to the comoving size 
(VArea) of about 6.7, 25.1 and 54.3 h^^Mpc at redshift 
0.1, 0.4 and 1, respectively, using the assumed cosmology. 
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The zCOSMOS is targeted to obtain - 20000 spectra of 
galaxies up to z < 1.4 with sampling expected to reach a 
rather uniform 60 — 70%. At the moment, we have avail- 
able spectra, and for most of them also possible various 
products, such as redshifts, luminosities, stellar masses, 
spectral line measurements etc. for 10644 objects, about 
half of the selected targets. This sample covers a slightly 
smaller area, of about 1.52 deg^. The spatial sampling of 
these objects is very inhomogeneous, with ~ 30% on av- 
erage. We refer to the catalogue containing only objects 
with reliable redshift measurements as the "10k sample" 
of the zCOSMOS survey. 

We build also the so called "30k sample" of galaxies in 
the zCOSMOS survey area which satisfies the magnitude 
selection criteria Iab < 22.5, but for which spectroscopic 
measurement is not yet available or is not of sufficient re- 
liability. For these galaxies we measure their photometric 
redshift probability fun ction P{z) using the ZEBRA code 
(jFeldmann et al.ll2006f ). The maximum in P{z) we sim- 
ply refer to as the photometric redshift. The uncertainty 
in the photometric reds hifts which we use is estima ted 
to be Sz = 0.023(1 -I- z) (jOesch et aLllin preparatioiH ) . 

The ZEBRA code has been employed also to calculate 
the rest-frame magnitudes of all Iab < 22.5 targets. The 
magnitudes are obtained as the best fit to the empirical 
set of spectral energy distribution (S EP) templates nor- 
malised to each galaxy photometry (jCapak et al.ll2007^ 
at the best available redshift (spectroscopic or photo- 
metric). The stellar masses are obtained by fitting stel- 
lar population synthesis models to the SED of t he ob- 
served magnitudes (,Bolzonella et al...2009;.Pozzetti et al.l 
[20091 ). 

2.2. Density field reconstruction 

For the various COSMOS projects, environment has 
been characterised in a few different ways. The full 
volume density field reconstruction h as b e en ca rried 
out by[Sco villc et al. (2007b), Masscv cCS] (|2007f l and 
iKovae et al.. (2009, ). .Scoville et al., (2007B) reconstruct 
the galaxy large scale structure at z < 1.1 using pho- 
tometric redshifts of galaxi es down to Iab < 2 5. Based 
on the observed shear field. iMassev et all (|2007f ) produce 
maps of the large scale distribution of dar k matter, re- 
solved in both angle and depth, up to z = 1. iKovac et al.l 
(|2009( l describe the reconstruction of the density field in 
the zCOSMOS volume up to z = 1, where the main 
ingredients are the high quality spectroscopic redshifts 
(with the uncertainty of about 100 km s~^) of galax- 
ies with Iab < 22.5 or sub samples of those galaxies. 
Moreover. IKovac et al.l (|2009( l present an overview of the 
all estimators of the continuous environment applied to 
the 10k zCOSMOS sample and discuss the importance 
of the scientific application on the exact choice of the 
density reconstruction method. In the following text, we 
will summarise the main steps in the zCOSMOS density 
field reconstruction method and define the galaxy sam- 
ples used for this process. We refer to Kovac ct al. (200i) 
for all the details. 

Briefly, the density p(r) at a given point in space 
r{RA, DEC, z), can be defined as: 



where mi is a weight related to some astrophysical prop- 
erty of an object (e.g. mass), W{\r — ri|;i?) is a spatial 
smoothing function, and 4){ri) is a function correcting 
the observed sample of objects to the complete sample of 
the same type of objects. The summation in Equation [8] 
goes over the observed galaxies which are used to recon- 
struct the density fi eld. We call them tracer galaxies. In 
IKovac et al.l (|2009| ) we discuss thoroughly the possible 
values/functional forms of m^, W{\r — ri|;i?) and (?!)(ri) 
and a choice of tracer galaxies and their influence on the 
reconstructed density field. 

We developed a new technique (ZADE) which allows 
us to incorporate both galaxies with spectroscopic and 
photometric redsh ifts in the density field reconstruction 
(jKovac et al.l2009| ) . This approach is based on modifying 
the P{z) function of galaxies from the 30k sample based 
on the proximity of galaxies from the 10k sample. We 
use counts of objects from the 10k sample in spheres of 
RzADE along redshift to redistribute the ZEBRA calcu- 
lated P{z) into a new, ZADE- modified, Pzade{z)- We 
have done extensive tests on mock catalogues to justify 
this method. For the zCOSMOS survey, a density field 
reconstruction with Rzade = 5 h~^Mpc produces a den- 
sity field without any major systematic error in the range 
0.1 < z < 1 (see the fo llowing Sect ion |43| and also Figure 
4 in IKovac et aLll2009f ). 

With the ZADE method, we use all Iab < 22.5 galax- 
ies in the zCOSMOS area to reconstruct the density field, 
albeit with different quality of their redshift measure- 
ment. When calculating densities of the tracer galaxies, 
tracers with spectroscopic redshift are counted as one 
(i.e. (5-function) at their measured redshifts, and the 
tracers with only photometric redshifts are counted as 
fractional objects according to Pzade(z) value at the 
redshift of consideration z (see also equation [Hj below). 
Formally, the ZADE approach is equivalent to 0(ri) — 1 
for every tracer galaxy. This means that in the ZADE 
approach the mean intergalaxy separation is a charac- 
teristic separation of the total population of that sample 
of tracer galaxies, and not only of a part of that pop- 
ulation with the measured spectroscopic redshift. The 
mean intergalaxy separation in the ZADE approach will 
therefore always be smaller than in some incomplete set 
of tracers, which is important in suppressing the shot 
noise effects. 

In our reconstruction method we assume that the 
structures traced by the 30k sample of tracer galaxies 
are delineated by the 10k sample of tracer galaxies (the 
sample with reliable spectroscopic redshifts), i.e. that 
the 30k and 10k samples are equally biased tracers of 
the underlying matter distribution. Based on the selec- 
tion criteria and observed distributions of the zCOSMOS 
galaxies, this assumption is valid up to z = 1 at least. We 
have done a number of careful checks on mock catalogues 
to validate this point (see Section |473|) . 

For the biasing analysis in this paper, the density of 
galaxies needs to match exactly the method with which 
the density of matter is estimated. Here, we count ob- 
jects within a spherical top-hat filter W{\r — ri|; R) with 
smoothing scale R = Rth , following 



(8) 
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W{\r-n\;R) 



d 



-Pi,ZADEiZi) if |r- Til <i?Tff 

otherwise, 
(9) 

In the equation above, |r — ri| is the three dimensional 
distance between a tracer galaxy and the point where the 
density is estimated. 

We express the reconstructed density field in terms of 
the dimensionless density contrast S{r) (commonly re- 
ferred to also as "overdensity" ) defined as 



S{r) = 



Pm{z) 



(10) 



Here Pm{z) is the mean density at redshift z, which we 
calculate as the volume average. For a point r in which 
a part of the cell defined by the smoothing scale R of 
a filter W falls outside of the survey limit, we apply an 
edge-correction by scaling the measured density with the 
fraction /e of the cell which is inside of the geometrical 
limits of the survey, such that the corrected density is 
pI fe- For the analysis, we use only the cells for which at 
least half of the volume is within the survey limits. 

The basic samples of all possible tracer galaxies (i.e. 
tracers selected only by Iab < 22.5), with both spec- 
troscopic and ph otometric redshifts , are equivalent to 
those described in lKovac et al.l ()2009l . Section 4.2.1). For 
the density field reconstruction we define four types of 
tracer galaxies: the "flux limited sample" of galaxies 
with Iab < 22.5 in 0.1 < z < 1 and three "luminos- 
ity complete samples" of galaxies satisfying the criteria 
Mb < -19.5 - z in 0.1 < z < 0.7, Mb < -20 - z in 
0.1 < z < 0.9 and Mb < -20.5 - z in 0.1 < z < 1. 
We include the passive evolution of AMb = — Az in se- 
lection of the luminosity complete samples in order to 
keep similar galaxies in a unique sample at every red- 
shift. The samples are selected to be complete for galax- 
ies of both red/ea rly and blue/late types (see Figure 1 in 
iKovac et al.ll2009[ ). 

We estimate pm{z) for each of these samples by adding 
up contributions of the tracer galaxies at each redshift. 
The contribution of any tracer galaxy to the mean den- 
sity at some z, up to Zmax, is calculated as ISV {z) /Vmax, 
where Ay(z) is the volume of the individual redshift bin 
and Vmax is the volume of the zCOSMOS survey in which 
a tracer galaxy of consideration can be detected, limited 
by Zmax- We take the passive evolution into account. 
The K-correction for the individual galaxies is obtained 
from the ZEBRA code. For practical purposes redshift 
is quantified in bins each Az — 0.002 wide. 

The redshift distribution of the mean volume densi- 
ties, weighted by unity, obtained by averaging the num- 
bers of zCOSMOS galaxies in 0.05 redshift bins and the 
corresponding smooth pm{z) functions are presented in 
Figure [TJ From this figure it is obvious that even the 
survey of the area of the zCOSMOS is dominated by in- 
homogeneities in the redshift distribution of galaxies of 
Az 0.1. On the other hand, the procedure to calcu- 
late Pm{z) described above smooths the peaky distribu- 
tion of the observed tracer galaxies, producing the mean 
volume density of luminosity complete tracer galaxies al- 
most constant with redshift. The mean intergalaxy sep- 
aration I of the zCOSMOS galaxies based on Pm{z) is 
~4.6, ~5.5 and -7.1 h^^Mpc for the Mb < -19.5 - z. 



Mb < —20 — z and Mb < —20.5 — z samples, respec- 
tively. For the flux limited sample of tracer galaxies, I 
increases from ~1.9 h~^Mpc at z = 0.1 to —3.3 h^^Mpc 
at z = 0.4 and to -6.7 h'^Mpc at z = 1. 

3. MEASURING THE BIASING FROM THE 
OBSERVATIONAL DATA 

A number of methods to measure the nonlin ear bias- 
ing fro m the observed data have been proposed. iSzapudil 
(Il998f ) suggests to use the cumulant correlators of the 
observed distribution of ga laxies in redshift surveys; 
iMatarrese et al.l ()1997D and IVerde et al] (|1998D derive 
the first two Taylor coefficients of the biasing function 
from the bispectrum of galaxies from the survey. SBD 
proposed a method, tested on simulations, to measure 
the mean biasing function relating the cumulative dis- 
tribution functions (CDFs) of the density fluctuations of 
galaxies and mass. 

In this work, we follow the method developed by SBD 
to derive the mean biasing function b{d) (or more specifi- 
cally the conditional mean function {6g\S)) and its second 
moments using the 10k zCOSMOS-bright survey over the 
redshift interval 0.4-1.0. Here, we summarise the main 
points and limitations of the SBD method and refer the 
interested readers to the original paper for all the details. 

We denote the CDFs of density fluctuations of galaxies 
and mass Cg{5g) and C((5), respectively. Both distribu- 
tions are functions of scale and redshift at which the 
density contrast is estimated and Cg{Sg) also depends 
on the type of tracer galaxies used to reconstruct the 
overdensity field. Assuming that the biasing relation be- 
tween galaxies and mass is deterministic and monotonic, 
the biasing function could be estimated from the inverse 
CDF of galaxies at a given value (percentile) of the CDF 
of mass (SBD): 



Sgi5) = C-'[C{S)], 



(11) 



where Cg ^[■■■] denotes the inverse CDF of galaxies. SBD 
used simulations to show that equation [TT] can success- 
fully reproduce the true conditional mean function {6g\6) 
over the broad range of S probed, overestimating it 
slightly in the most dense regions. They also conclude 
that the assumption of the monotonicity of {Sg\S) is valid 
(or represents a very good approximation) on scales of a 
few Mpc, which are the scales over which the biasing is 
measured in practice. 

3.1. PDF of density contrast of galaxies and mass 

Current galaxy redshift surveys provide data to cal- 
culate the galaxy density contrast and its PDF with 
suflicient accuracy u p to redshifts of about 1.5 (e.g. 
iMarinoni et al.l [20051 ). Obtaining the mass density con- 
trast observationally is on the other hand an extremely 
difficult task. Direct reconstruction has been possible 
using galaxy peculiar velocities as tracers of mass fiuc- 
tuations so far only in the l ocal universe (out to — 100 
h-iMpc; iDekel et al.l 1199(1. [l999l: IPek el l200l or using 
weak lensing up to z — 1 (jMassev et a l. 2007), although 
with much lower resolution in the redshift dimension 
than is obtained for the galaxy density field from the 
galaxy redshift surveys. From a theoretical point of view 
it has been shown that in comoving space, the density 
contrast S of matter follows, to a good approximation. 
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a lognormal distribution p{6) (| Coles fc Jonesl[l991[ ). In 
this paper we will use this theoretical approximation of 
the matter distribution and in the following text, we will 
summarise the theoretical development necessary to cal- 
culate the PDF of the matter. However, we want to 
mention that there is a project in development to use si- 
multaneously galaxy and matter density field (the latter 
one reconstructed from the weak lensing shear maps) in 
the COSMOS volume to estimate bias directly from the 
observed data. 

The lognormal distribution of the matter fluctuations 
can be expressed as 



[ln(l-f (5)+w|/2]2 



(12) 



In the last equation, the parameter uj'^ is defined as 

ujl = \n[l + {5^)R]. (13) 



where 5 is directly related to the variance <jj^{z) of the 



density contrast field at redshift z via 



(14) 



given that the density contrast field has zero mean. In 
the equations above, the index R denotes the smooth- 
ing scale at which the density field is reconstructed. The 
value of ub. on a given scale is determined by the adopted 
cosmology. For the smoothing scales R which are large 
enough to be in the linear regime, its evolution with red- 
shift can be modelled as 



or{z) = Gn{z = O)i^(z) 



(15) 



where D{z) is the linear growth factor of density fluctu- 
ations, normalised to unity at z = 0. For smaller scales, 
one would need nonlinear corrections. 

Given that the derived lognormal form of the PDF of 
the matter overdensity field is calculated in real (comov- 
ing) space, and that the PDFs from surveys are obtained 
in redshift space, one has to convert both functions to 
the same space. If the redshift distortion affects both 
galaxies and matter in a similar way, we can expect that 
the mean biasing function 6(5) and the conditional mean 
function in z space will be similar to the one in 

real space: 



(16) 



SBD show that the PDF (or CDF) of the mass density 
contrast in redshift space can also be well described with 
the lognormal function (equation [12]), with standard de- 
viation obtained in redshift space. The relation between 
the azimuthally averaged standard deviations of mass 
fluctuations in real (7r{z) and redshift crf^(z) comoving 
space at high redshift is (e.g. SBD) 

aUz)-[l + lfiz) + lfiz)]'/'an{z) (17) 



based on the expression derived by iKaiserl (|1987) . For 
the calculations, we use the relation between the growth 
rate / and the growth factor D given by 



f{nm{z),nA{z)) 



dlnD 



+ f7A(z) + 



rflna 



= -1 - 

2 

where a — {1 + z)^^ and 

D(z) = , 

' gmi + z) 

for linear fluctuations (jCarroll et al.lll99^ and 



(18) 



(19) 



g{z) « -r2„(z) X 



nt/'iz)^nA{z) 



1 



1 



»a(z) 
70 



' (20) 



In the literature, there are other approximations avail- 
able for the growth rate / at different redshifts. For 
example. iLahav et al.l (IT99ll ) derived / « 17™ (z)"-^ and 
I Wang fc Steinhardtl (11998!) derived / « n^izf-^^, the 
latter one allowing the possibility that the current ac- 
celerated phase of the universe is due to quintessence (a 
time- evolving, spatially inhomogeneous component with 
negative pressure) and not due to the cosmological con- 
stant. 

4. ERRORS IN THE MEAN BIASING FUNCTION AND 
BIASING PARAMETERS 

The mean biasing function (or the conditional mean 
function {Sg\6)) and its moments derived using the 
method described in Section [3] will contain uncertainties 
due to the finite volume of the survey (cosmic variance 
errors) and the use of discrete tracers (i.e. galaxies) to 
reconstruct the continuous density field (shot noise er- 
rors). In addition to this, the results on the biasing will 
contain errors related to the reconstruction of the density 
field. We use mock catalogues to asses the contribution 
of each of these errors. 

The mock catalogues which we employ are based 
on the light cones for the COSMOS survey kindly 
provided by iKitzbichler fc White! (|2007i) . The mock 
catalogues a re built from the N-body Millennium 
Simulation (|Springel et al.l !2005! ) using the semi- 
analytic modelling of g alaxy properties based on 
!De Lucia fc BlaizotI (!20071) wit h some modifications de- 
scribed in !Kitzbichler fc White! (|2007! ). We cut the mock 
catalogues such that their area and redshift limits match 
the zCOSMOS survey. Depending on the errors which 
we test, we use mocks with different flux or luminosity 
limits and of different sampling. 

We reconstruct the overdensity field on a grid with 
separation of 0.5 h~^Mpc in the RA — DEC plane and 
0.002 in the z direction. As for the real data, we use ev- 
ery tracer galaxy to reconstruct the density field, using 
either its spectroscopic redshift or its ZADE-modified ar- 
tificial photometric redshift probability distribution. As 
we have mentioned previously, this is equivalent to (f>i = 1 
for every tracer galaxy following notation of equation |8l 
We limit ourselves to the use of the unity-weighted over- 
densities, i.e. TTii = 1 in the same equation. The mean 
density in the individual mock catalogues is estimated 
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by smoothing the N{z) distributions for the flux hm- 
ited samples. For the luminosity complete samples the 
mean density is estimated for each redshift bin simply 
as the number of objects divided by the volume. While 
for the biasing analysis we use redshift bins of 0.3 width, 
to estimate the mean number of objects we use redshift 
bins of 0.4. These bins are broadened symmetrically by 
Sz = 0.05. 

To be consistent with the analysis of the zCOSMOS 
data presented in Section [51 when calculating the mean 
biasing function using the mocks (in order to evaluate 
various errors) , we assume cs = 0.8, roughly consis - 
tent with the latest WMAP results (jPunklev et al.ll2009f ). 
even though the mock catalogues are based on the simu- 
lations with (78 = 0.9. Similar to recent works, we choose 
the logarithmic representation of overdensities in order to 
emphasise the behaviour in the underdense regions. For 
simplicity, we will refer to the conditional mean function 
of log(l + Sg) at a given log(l + 6) also as {Sg\S). 

Before we proceed with measuring the errors in the 
conditional mean function {Sg\5) function and the second 
moments of the mean biasing function, we need to esti- 
mate the scale R (smoothing top-hat window) at which 
we can reconstruct the overdensity field without serious 
systematic errors. We reconstruct the galaxy overdensity 
field on the so called 10k-|-30kZADE and 40k mock cat- 
alogues. Both of these types of mock catalogue have the 
same geometrical constrains as the zCOSMOS survey. 
In the 40k mock catalogue all galaxies have measured 
spectroscopic redshifts. In the lOk-f SOkZADE catalogue, 
only ^ 10000 galaxies with Iab < 22.5 have known spec- 
troscopic redshifts. Their RA — DEC distribution has 
be en matched to th e complicated sampling (see Figure 4 
in lLillv et al.ll2069f l and the redshift success rate of the 
10k zCOSMOS sample. We take into account Iab < 22.5 
objects without redshift by defining their photometric 
redshift probability distribution P{z). In the mock cata- 
logues, we model the initial P{z) to be a Gaussian distri- 
bution with (T = 0.023(1 -|- z), randomly offset in redshift 
using offsets selected from the same distribution. For the 
density field reconstruction this probability is then mod- 
ified using the ZADE algorithm to yield Pzade{z) (see 
Subsection [221) ■ For computational purposes, PzADsiz) 
is discretised in redshift bins of 0.002. 

To test the minimum smoothing scale, we use Iab < 
22.5 samples of galaxies and three smoothing scales: 5, 8 
and 10 h^^Mpc. As an example, the PDFs of the recon- 
structed density contrast for the 40k and 10k-t-30kZADE 
samples is presented in Figure [2l for one of the mock cat- 
alogues. For the current status of the zCOSMOS survey 
we need scales of at least 8 h~^Mpc to reconstruct the 
overdensity field with acceptable errors at every Sg up to 
z~ 1. 

4.1. Cosmic variance 

One of the errors entering the biasing analysis occurs 
due to the relatively small survey volume and consequent 
noise in Cg{dg). This error is usually termed cosmic vari- 
ance. To quantify it, we employ 12 mock catalogues de- 
signed to be equivalent in the terms of geometrical and 
magnitude selection constrains to the zCOSMOS survey. 
We use the full 40k catalogues to select the different sub- 
samples of tracer galaxies to reconstruct the density field. 

As an example, we show in Figure [3l the conditional 



mean functions {Sg\S} obtained for the 12 mock cata- 
logues, where the mock galaxy density field has been re- 
constructed with the fiux limited sample of galaxies and 
top-hat fiher of 10 h^^Mpc in OA < z < 0.7. Due to 
cosmic variance, there is a dispersion of Sg values corre- 
sponding to every S value in the biasing analysis since 
every mock will uniquely map 5 to a single (mean) 6g 
value. In the mocks, the measured range of possible Sg 
values is largest in the most underdense regions, where 
the standard deviation a of log(l + Sg) values is a ^ 0.1. 
The (5g-dispersion gets smaller towards regions near the 
mean density, where a < 0.05. In the most overdense 
regions, the cosmic variance error is increasing again, up 
to a ^ 0.05. We will use the standard deviation in the 
conditional mean function {Sg\S) and in the biasing pa- 
rameters obtained from the 40k mock catalogues as the 
errors in corresponding values derived from the zCOS- 
MOS sample. 

It is important to point out that the the conditional 
mean function {Sg\S) in each of the mock catalogues is 
derived with respect to the same theoretical underlying 
matter PDF, as it would be the case with the observa- 
tions. To obtain the "true" cosmic variance, one would 
need to derive {Sg\S) with respect to the PDF of the 
matter using the DM particles in each of the mock cat- 
alogues (which we do not have). Moreover, stochasticity 
may also contribute to the measured scatter in the {Sg\S) 
function. 

4.2. Shot noise errors 

Sparse sampling artificially enhances both the positive 
and the negative tails of the density contrast distribution. 
This effect broadens the PDF and thus steepens the con- 
ditional mean function {Sg\S) (e.g. SBD). To evaluate 
the shot noise errors, we first use the 12 mocks including 
the full geometrical constrains of the zCOSMOS survey 
and limited by the rest-frame magnitude Mb < — 18 (ig- 
noring the luminosity evolution for the moment). This is 
the faintest magnitude cut for which the given mock sur- 
veys r < 26 (jKitzbichler fc White! 120071 ) are luminosity 
complete up to z = 1. The mean separation / between 
galaxies is about 2.7 h~^Mpc in the individual mock cat- 
alogues, which should be sufficient to obtain the mean bi- 
asing function and its moments without shot noise errors 
for scales of 8 h~^Mpc and larger. To quantify the effect 
of the shot noise errors we then resample these mock 
catalogues at random such that the mean distance be- 
tween galaxies I is 3, 4, 5, 6, 7 and 8 h^^Mpc. We derive 
the conditional mean function {Sg\S) and the second mo- 
ments of the mean biasing function using the overdensity 
fields reconstructed from the resampled mock catalogues 
in two redshift bins: Az = 0.4 - 0.7 and Az = 0.7 - 1. 

The effect of the sampling on the {Sg\S) function de- 
rived from the zCOSMOS size survey is presented in Fig- 
ure m for one of the mock catalogues. The shot noise 
errors can change the shape of {Sg\6) in the underdense 
regions dramatically, the effect clearly depending on the 
value of I (relative to i?), i.e. on then mean number of 
galaxies per sampling region. In the underdense regions, 
the value of the mass density contrast associated with 
a given galaxy density contrast is artificially shifted to 
higher values when subsequently sampling smaller frac- 
tions of the population of tracer galaxies. This is im- 
portant given that this minimum S value, below which 
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the galaxy density field does not trace that of the mat- 
ter, is often interpreted as the minimum mass density 
contrast below which the formation of tracer galaxies is 
partially or completely s uppressed fe.g. ISomerville et al.1 
l2001l : lMarinoni et al.l f2005). 

The summary of the effect of the shot noise errors 
on the second moments of the mean biasing function is 
shown in Figure^ We average the parameters b and b/b 
estimated from the 12 mocks and present the averaged 
parameters with the errors calculated as the rms values 
of the corresponding parameters in these 12 mocks. The 
linear biasing parameter b increases constantly with in- 
creasing mean intergalaxy separation, while there is some 
indication that the nonlinearity parameter b/b slightly 
decreases with increasing mean intergalaxy separation, 
at least in the higher redshift bin. For example, for the 
smoothing scale 8 h^^Mpc the linear biasing parameter 
increases by ~ 20% and nonlinearity decreases by ~ 0.6% 
when changing the mean intergalaxy separation from 3 
to 8 h~^Mpc. The effect of the shot noise errors on the 
nonlinearity parameter can be practically neglected on 
smoothing scales of 8 h^^Mpc and larger. Note that the 
presented standard deviations of the biasing moments 
from different mocks are dominated by the cosmic vari- 
ance. 

4.3. Reconstruction errors 

The reconstructed galaxy density field is meant to be 
the density field of the full population of tracer galax- 
ies. In practice, the whole population of galaxies is 
never available, and the measured properties of sample 
galaxies, relevant for the reconstruction, contain some er- 
rors, such as errors in spectroscopic redshifts or ZADE- 
modified redshift photometric probability functions in 
the samples of tracer galaxies. We refer to reconstruc- 
tion errors as those errors which arise due to the fact that 
only a part of the total population of tracer galaxies with 
their measured properties is used for the density field re- 
construction instead of the full population with their true 
properties. Given that we use the ZADE approach, we 
can easily differentiate between the reconstruction and 
shot noise errors, as defined here. However, when using 
only a fraction of tracer galaxies with measured redshifts 
to reconstruct the density field (and correcting statisti- 
cally for the galaxies without measured redshifts), it can 
be difficult to separate between such defined reconstruc- 
tion and shot noise errors. 

To gain an understanding of the effect of reconstruction 
errors in the 10k zCOSMOS density field on the biasing 
analysis, we derive the conditional mean function {Sg\S) 
and biasing parameters using 12 mock 10k-|-30kZADE 
catalogues, for the two samples of tracer galaxies: flux 
limited Iab < 22.5 and volume limited Mb < — 20 — 
z. We compare the results to the corresponding ones 
derived from the 40k mock catalogues. We use only the 
smoothing scale of 8 h^^Mpc, given that the error in the 
density field reconstruction is gene rally larger for sma ller 
smoothing scales (see Section 5 in lKovac et al.l[2009[ ). 

The resulting {Sg\S) functions for the lOk-^SOkZADE 
(observed) and 40k (true) mocks for smoothing top-hat 
scale of 8 h^^Mpc, obtained by averaging the conditional 
mean functions {Sg\S) from the individual mock cata- 
logues, are presented in the lower panels of Figures [6] 



and [7] for the flux limited and Mb < —20 — z limited 
samples of galaxies, respectively. The error in the re- 
construction of the density field changes the shape of 
{Sg\S) in different ways for the flux and luminosity com- 
plete samples of tracer galaxies, as following. The re- 
construction errors for the flux limited sample of galax- 
ies are negligible in low redshift bins, and they increase 
in higher redshift bins, artificially increasing local bias 
b{S,z,R) values in the regions of the both highest and 
lowest density contrasts. For the volume limited samples 
the reconstruction error is manifested by lowering the lo- 
cal bias b{S, z, R) values in the most underdense regions, 
and this error decreases with redshift. The different ef- 
fect of the reconstruction error in the fiux and luminosity 
complete samples arises from the fact that we are using 
the same ZADE-modified probability functions for the 
objects without spectroscopic redshift in both samples, 
but the overall numbers of galaxies in the flux and lumi- 
nosity complete samples are different. These differences 
change with redshift (see Figure [T]) , and therefore the 
reconstruction error in the biasing analysis manifests it- 
self differently in different redshift bins. In summary, 
the reconstruction error is most notable in the S < 
regions. At the reconstructed value of galaxy density 
field log{l + 5g) — —1, the reconstruction error can cause 
differences up to ~ log(l + 6) = 0.1. Any change in 
the conditional mean function {Sg\S) in the underdense 
region smaller than this value can therefore not be re- 
garded as significant. 

In the upper panels of Figures [6] and [7] we present the 
corresponding rms values (a) of the average {6g\5) value 
from these 12 mocks, over the whole range of S values. 
These a values correspond to the cosmic variance in the 
reconstructions of the galaxy density field with the 40k 
and 10k-f-30kZADE samples. For log(l + S) > -0.5, 
a is almost identical for the 40k and 10k-|-30kZADE 
reconstructions for the both fiux and luminosity com- 
plete samples. The fact that the conditional mean func- 
tion {Sg\S) does not vanish at log(l + 5) < —0.5 for 
the lOk-f 30kZADE reconstruction with the luminosity 
complete samples is due to the ZADE-approach. Given 
that the matter distribution is lognormal, the total effect 
of this artificial "filling" of the underdense regions with 
"fractional" galaxies is negligible on the final results. 

In addition, we compare the b and b/b parameters from 
the 10k+30kZADE and 40k samples (Figure (5]). The re- 
construction error in the b parameter acts as an artificial 
increase of this parameter for both fiux limited and lu- 
minosity complete tracer galaxies. On the other hand, 
we see that the nonlinearity parameter b/b is almost not 
affected by the error in the reconstruction. 

5. BIASING ANALYSIS: RESULTS AND THEIR 
INTERPRETATION 

We derive the conditional mean function {Sg\S) and the 

biasing parameters b and b/b of the 10k zCOSMOS galax- 
ies using four galaxy samples to reconstruct the den- 
sity field and obtain the PDF of the density contrast: 
Iab < 22.5 fiux limited sample and three luminosity 
complete samples of Mb < —19.5 — z, Mb < —20 — z 
and Mb < —20.5 — z galaxies. We derive the density 
field following equation [S] with the top-hat three dimen- 
sional filter (equation [9]) , using smoothing scales of 8, 
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10 and 12 h^^Mpc. The range of the smoothing scales 
is Umited by the minimum of 8 h~^Mpc at which we 
can rehably reconstruct the overdensity field (based on 
the reconstruction method tested on the mocks) and the 
transverse size of the zCOSMOS field. We use the ZADE 
approach to account for galaxies without reliably mea- 
sured spectroscopic redshift and therefore = 1 for ev- 
ery tracer galaxy. If not otherwise stated, we use the 
unity- weighted overdensity field (mj = 1) for the biasing 
analysis. The zCOSMOS selection catalogue is based on 
the 0.1 arcsec resolutio n HST images in F814W filter 
(|Koekemoer et al.ll2007[ ). supplemented by photometry 
from a high resolution CFHT image in i filter in the case 
when the HST data is missing or when the HST photom- 
etry is not of required quality. With this combination of 
underlying imaging only a negligible part of the zCOS- 
MOS area is affected by the foreground stars (|Lillv et al.l 
l2009f ). and there is no need for masking of some regions 
for the density field reconstruction. 

As for the reconstruction of galaxy overdensity field 
on the mocks that was used to obtain the contribution 
of various errors on the conditional mean function((5g 
and the second moments of the mean biasing func- 
tion (described above), we reconstruct the zCOSMOS 
galaxy overdensity field on a grid with the points in the 
RA — DEC plane separated by 0.5 h~^Mpc and with 
Az = 0.002. The overdensity field reconstructed with 
R — 8 h~-'^Mpc and flux limited tracers in 0.4 < 2 < 1 is 
presented in Figure [31 The complex, cosmic- web appear- 
ance of the density field, consisting of cluster-like struc- 
tures, surrounded by empty, void-like r egions, is detected 
in the whole redshift range probed (see lKova6 "eran[2009l 
for the more detailed discussion of the structures in the 
zCOSMOS overdensity field). 

We carry out the biasing analysis in the redshift range 
0.4 < z < 1, starting at z = 0.4 in order to exclude 
redshift slices in which for the majority of grid points the 
large fraction of filter aperture (i.e. larger than 0.5) falls 
outside of the survey volume. We calculate the density 
contrast in four redshift bins: 0.4 < z < 0.7, 0.5 < z < 
0.8, 0.6 < z < 0.9 and 0.7 < z < 1, which overlap in 
order to suppress the effect of the cosmic variance. This 
is clearly visible in the zCOSMOS ove rdensity field plots 
(FigurelSl see also Figures 12 and 17 in lKovac et al.ll2009( ) 
and hampers the expected evolution in the distribution 
of lar ge scale structures with cosmic time (see Figure 
19 in iKovac et"al] l2009f ) . At lower redshift, the mean 
overdensity in the individual redshift slices differ from 
zero by a few percent. However, in 0.6 < z < 0.9 and 
0.7 < z < 1, the mean density is about 0.1 — 0.15 for all 
but Mb < —20.5 — z sample in 0.6 < z < 0.9, for which 
the mean overdensity is about 0.2. Therefore, the bias 
values for these samples needs to be taken with caution. 

As discussed in previous sections, the PDF of the mass 
density contrast is calculated assuming a lognormal dis- 
tribution (equation [12]), specified by the adopted cosmol- 
ogy. The values of measured biasing parameters are pre- 
sented in Table [TJ The detailed discussion and interpre- 
tation of the results obtained from the biasing analysis is 
presented in the following subsections. When interpret- 
ing the zCOSMOS biasing results and comparing them 
to other work, one has to keep in mind that the obtained 
results and their exact redshift evolution is cosmology 
dependent. This is evident from the dependence of the 



growth rate on the cosmological parameters, and partic- 
ularly the results are dependent on the chosen as{z — 0) 
normalisation. SBD find that the linear biasing param- 
eter b changes as cr~^, while the nonlinearity parameter 
b/b changes only very weakly with a, b/b ^ a'^-^^. We 
show in Figure [10] the change of the shape of the mean 
biasing function with the ag parameter. Lowering erg 
will produce less structure at a given epoch, and it has 
a similar effect on the change in the shape of the mean 
biasing function as the increasing the mean intergalaxy 
separation. 

5.1. Shape of the conditional mean function {Sg\S) 

We show the conditional mean function {Sg\S) of galax- 
ies in Figures [Tl] and 1121 where the galaxy overdensity 
fields are reconstructed with the luminosity complete 
Mb < —20 — z sample of galaxies for the smoothing 
filters of 8 and 10 h^^Mpc, respectively. We also present 
the corresponding linear biasing approximation Sg = b^S 

with b^ — b at every S. 

The {6g\6) function vanishes in the most underdense 
regions. At moderate underdensities the {Sg\S) function 



sharply rises as we approach Sg 



0. Starting in 



these and up to the mildly overdense regions, the {6g\6) 
function closely follows the linear relation with 6. In the 
most overdense regions, our results suggest that galaxies 
are antibiased tracers of the underlying matter distribu- 
tion. The local slope b{5, z, R) of the biasing relation in 
the underdense regions is larger than unity. In the over- 
dense regions, the trend is less clear, as the local slope 
can take values both larger and smaller than unity. 

This characteristic shape of the conditional mean func- 
tion {Sg\S) persists for all the samples of tracer galaxies 
and redshift intervals covered by our study. The bias- 
ing relation between galaxies and matter is clearly non- 
linear in the most underdense and overdense regions in 
0.4 < z < 1, in agreement with previous work at these 
redshifts based both on simulations and semi- analytical 
modelling (SA M; SBD, ISomerville et al.ll200i[ ) or obser- 
vations (Mari noni et al.l 2005( 1. 

Theoretical works provide some explanation for the 
observed shape of the conditional mean function {Sg\6). 
Vanishing of the function in the underdense regions can 
be interpreted within a scenario in which galaxies do not 
form below some mass density threshold. However, as 
discussed above, one needs to be careful in this inter- 
pretation and take the possible shot noise effects into 
account before inferring this threshold. The antibiasing 
of galaxies in the most overdense regions observed at red- 
shifts of unity and below can be explained by quenching 
of galaxy formation in these regions, as at these redshi fts 
the densest regions become too hot (|Blanton et al.|[200QV 
The other possibility to explain the antibiasing of galax- 
ies in the most overdense regions are different epochs of 
forma tion of galaxies in ove rdense and underdense re- 
gions (jYoshikawa et al.ll2001f ). This follows fr om the hi- 
erarchical scenario of galaxy formation, where for a given 
mass scale, there is a tendency for objects in overdense 
regions to form earlier than objects in the underdense re- 
gions. Therefore the young galaxies (e.g. w ith formation 
redshift since 1.7 in Yoshika wa et ani2001f ) are expected 
to form in low density, which are also low temperature, 
environments. Also, the merging of galaxies in high 
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density environments could lower the nu mber density of 
galax ies used to derive the density field (|Marinoni et al.l 
[2005h . 

5.2. Scale dependence 

To examine the possible dependence of the shape of 
the conditional mean function {5g\6) on the smoothing 
scale, we present in Figure [13] {Sg\S) derived with the 
Mb < —20 — z tracer galaxies smoothed with three dif- 
ferent filters: 8, 10 and 12 h-^Mpc. The {Sg\S) func- 
tions for different smoothing scales are almost identi- 
cal to each other for the largest part of the explored 
range of overdensities, in both redshift bins which we 
plot (0.4 < z < 0.7 and 0.6 < z < 0.9). In the most 
overdense regions the characteristic S at which galaxies 
become antibiased shifts from higher to lower values as 
the smoothing scales shift from lower to higher values. 

We summarise the impact of the different smoothing 
scales on the biasing parameters in Figure [HI using the 
same overdensity field. We do not detect any signifi- 
cant dependence of the derived linear biasing parameter 
nor the nonlinearity parameter on the smoothing scale. 
However, we cover only a very narrow range in scales 
and taking the estimated errors into account, we would 
detect a dependence of the bias on the smoothing scale 
only if the effect is very strong. 

The negligible dependence of the biasing function and 
parameters on the scales of 8 h~^Mpc and larger is in 
agreement with the arguments from a number of the- 
oretical works, according to which bias is expected to 
be consta nt on scales l arger than a few h~^Mpc (e.g. 
Kauffma nn et all Il997t iMann et al] Il998t iBenson et all 
200Q). How e ver, using hydrodynamical simulations 
Blanton et al.l ()1999i^ detect decreasing of galaxy bias 
with the smoothing scale, where bias is defined as the 
ratio between the variances of the number of galaxies 
and mass withing spheres of radius R. This depen- 
dence is significant on the scales R smaller than the 
transition scale betwee n the linear and nonh near regimes 
(and it is 16 h ~ ^Mpc in lBlanton et al.lll999l simulations). 
iBlanton et al.l (|1999f) explain the scale dependence of the 
bias to follow from the dependence of the galaxy density 
field on the local temperature, which reflects the gravita- 
tional potential related to the mass density field. If the 
gas is too hot, galaxies will not form, infiuencing directly 
the reconstruction of the galaxy density field. 

5.3. Luminosity dependence 

In this subsection we investigate the dependence of the 
{Sg\6) function and biasing parameters on the luminos- 
ity of galaxies used to reconstruct the overdensity field, 
for a fixed smoothing scale of 10 h~^Mpc. The {dg\S) 
function derived with the tracer galaxies of different lu- 
minosity thresholds in two redshifts bins (0.4 < z < 0.7 
and 0.6 < z < 0.9) is presented in Figure [TCI We detect 
a weak dependence of {Sg\S) on the luminosity of the 
tracer galaxies. This effect is clearest in the overdense 
and mildly underdense regions of the mass density con- 
trast distribution, where the local bias of more luminous 
galaxies is higher than the local bias of less luminous 
galaxies. In the most underdense regions the differen- 
tiation between the conditional mean functions {Sg\6) of 
the galaxies with different luminosity thresholds becomes 



barely visible; however, in this regime {6g\S) of less lu- 
minous galaxies is systematically above {6g\S) of more 
luminous galaxies, although hard to see. Considering 
the reconstruction errors, the differentiation between the 
{Sg\S) functions due to the luminosity of tracer galaxies 
should be even more pronounced. We have seen in Fig- 
ures [Hi and [7| that due to the reconstruction errors the 
conditional mean function {dg\5) in the regions of the 
lowest mass density contrast is getting shifted in the op- 
posite directions for the fiux limited and luminosity com- 
plete samples, making galaxies to appear higher or lower 
biased for these two samples, respectively. Taking these 
errors into account, we expect the difference between the 
critical S below which Sg of galaxies do not trace the mat- 
ter any more to be about log(l -f (5) ~ 0.1 for the {5g\S) 
function of the flux limited and Mb < —20 — z samples 
of galaxies. 

The {Sg\S) function in lMarinoni et all (|2005| ) shows the 
same behaviour with luminosity in the most overdense 
regions, even though their errors in the most overdense 
regions are larger. Interestingly, iMarinoni et al.l (|2005f ) 
detect a much clearer differentiation of the conditional 
mean function {6g\S) with luminosity in the regions of 
(5 < 0, in the same direction as our data indicate. They 
interpret the change in the shape of the {Sg\S) function in 
the underdense regions with luminosity as an indication 
that the formation efficiency of galaxies is shifting to- 
ward s higher densities with increasing luminosity. How- 
ever, [MarinonreF^n (|2005f ) do not take into account the 
change in the mean intergalaxy separation between dif- 
ferent populations of galaxies, as this fact by itself will 
cause this type of shift in the {Sg\S) function (due to the 
different shot noise manifestation for the different popu- 
lations of galaxies as shown in Section |4?2|) . 

The mean biasing parameter measured by 6, presented 
in the upper panel of Figure [TH shows a systematic de- 
pendence on the luminosity of the tracer galaxies: it is 
higher for the more luminous galaxies. The errors in the 
plot include only the cosmic variance and they indicate 
the expected spread in the obtained b values at a given 
redshift. The errors from the reconstruction can be ne- 
glected at lower redshifts, but for z > 0.7 they artificially 
increase the b parameter for both the flux limited and lu- 
minosity complete samples (for about ^ 5% and ~ 4% 
at z ~ 0.75 for the flux limited and Mb < —20 — z 
luminosity complete samples, respectively, using a top- 
hat smoothing of 8 h~^Mpc, see Figure [8]). In addition, 
more luminous samples of galaxies will have larger mean 
intergalaxy separation, what will artificially increase the 
b parameter. For example, as seen in the mocks, the rela- 
tive increase in b (with respect to the results with Z ~ 2.7 
h-^Mpc) would be ~ 4.5% and ~ 10.5% at z ~ 0.75 for 
AIb < —20 — z and Mb < —20.5 — z samples of galaxies, 
purely because the differences in the mean intergalaxy 
separations. Given that we measure larger differences 
between the b parameters from different luminosity com- 
plete samples than expected by purely taking into ac- 
count mean intergalaxy separations for these different 
populations, we are confident that the trend of increas- 
ing b with luminosity refiects intrinsic physical processes 
in galaxy formation. The dependence of the correlation 
function on the luminosity is a known observational re- 
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suit (e.g. IZehavi et al]|2005t iPollo et al.ll2006l : [Coil et all 
l2008t but see lMeneux et al.ll2009l ). therefore an increase 
in the linear bias with luminosity is expected. 

We do not detect any significant dependence of nonlin- 
earity on the luminosity of the tracer galaxies. Based on 
the mocks, the errors from the density field reconstruc- 
tion can be neglected in the nonlinearity parameter. The 
shot noise effects in the nonlinearity parameter tend to 
be more important at higher redshifts. For example, the 
shot noise effect would be responsible for about 0.3% de- 
crease in the nonlinearity for the mean intergalaxy sepa- 
ration 7-8 h^^Mpc. 

5.4. Redshift evolution 

The evolution of the conditional mean function {Sg\S) 
and biasing parameters with redshift can be studied in 
all results presented so far. At a given luminosity and 
smoothing scale, the basic shape of the {Sg\S) function 
(e.g. Figures [TT] and [T2|) is preserved in all three red- 
shifts bins probed by the overdensity field traced by the 
10k zCOSMOS galaxies. The characteristic value of 5 
at which galaxies (apparently) do not trace matter any 
more, defined here by Sg < —0.9, shifts from S « —0.70 at 
z ~ 0.55 to (5 « -0.63 at z - 0.75 for the Mb <-20-z 
tracers. However, this difference is comparable to the re- 
construction error in the {Sg\S) function with this sample 
(Figure [7]), therefore we can not distinguish whether the 
redshift evolution in this characteristic 6 value is the real 
effect or not. On the other hand, at positive 5 galaxies 
become antibiased at gradually higher matter overdensi- 
ties when going from lower to higher redshifts. This value 
of 6 shifts from 4.46 at z ~ 0.55 to 6.75 at z - 0.75. 
Given that we have selected partially overlapping red- 
shift slices in order to reduce cosmic variance, we have 
also smoothed the possible redshift evolution in the shape 
of the {dg\S) function. 

The redshift evolution of the mean biasing parameter 
b is more evident (Figures [T3] and [TB|) . For the adopted 
cosmology, b increases from 1.24 ±0.11 at z ^ 0.55 to 
1.63 ± 0.15 at z - 0.75 for the Mb < -20 - z sample of 
galaxies for the top-hat smoothing of 10 h~^Mpc, taking 
very similar values for the other smoothing scales. Also, 
b is higher for the more luminous samples of galaxies (but 
see the discussion on errors in the previous Subsection). 

The nonlinearity of the mean biasing function as mea- 
sured by b/b is the least affected by the shot noise and 
reconstruction errors. The b/b do not show any signifi- 
cant redshift evolution, in addition to no dependence on 
the smoothing scale or on the luminosity threshold of 
galaxies used to reconstruct the overdensity field. The 
nonlinearity parameter is of the order of a few percent 
(~ 2%). 

The redshift evolution of the shape of the {dg\S) func- 
tion in the manner indicated by our results has bee n seen 
in both the simulations and SAMs (e.g. j S^omcrvill e et al.l 
[200 ll) and observations (M arinoni et al.f l2005). The red- 
shift evolution of the biasing is visible most clearly using 
the linear biasing parameter. The increase in the linear 
biasing parameter with redshift seen in the zCOSMOS 
data supp orts the the o retica l predictions of the bias evo- 
lution. Blanto n et al.l (|2000f ) discriminate between three 
different effects which are responsible for this bias be- 
haviour. First, galaxy formation process shifts from the 



highest peaks in the density field at early epoch to the 
lower peaks in the density field as time progresses, i.e. 
this process shifts from the most biased to less biased 
tracers of the underlying density field. Second, the for- 
mation of galaxies in the most dense environments is 
halted towards lower redshifts, because these regions be- 
come filled with gas which is shock heated and virialised, 
but which is not able to cool and collapse. At higher 
redshifts {z ~ 3) galaxies are expected to be biased even 
in the highest density regions, because they are still suffi- 
ciently cold enough to allow for fast cooling of the clumps 
of gas. These two effects combined are responsible for the 
shift in the galaxy formation from the most dense to less 
dense environments. Third, once galaxies are formed, 
they experience the same gravitational physics as the 
dark matter, and therefore the distribution of galaxies 
and matter becomes mor e and more similar ()Frvlll996t 
iTegmark fc Peeble£lll998[ ). 

5.5. Biasing for the luminosity or stellar mass weighted 
density field 

The general equation to reconstruct the galaxy den- 
sity field. Equation [51 is commonly used by weighting all 
galaxies in the tracer sample equally, i.e. using TOj = 1. 
However, it can be expected that some galaxy properties, 
particularly total mass, are better tracers of the underly- 
ing matter density field. To address this quantitatively, 
we repeat the galaxy density field reconstruction in the 
zCOSMOS volume using the sample of Mb < —20 — z 
and AIb < —20.5 — 2 tracer galaxies as defined in Sub- 
section I2.2[ but now using two different mi-weighting 
schemes. We use = Lba and rrii = M*_i, which 
are estimated reasonabl y well for practically all ga lax- 
ies in the 40k sample fsee lOesch et al.|[in preparation! and 
iPozzetti et al.|[2009l for the absolute magnitude and stel- 
lar mass measurements). 

The conditional mean function {6g\S) between the Lb- 
and M:,, -weighted galaxy overdensity field and the mat- 
ter overdensity field for the AIb < —20 — 2 sample of 
tracer galaxies is presented in Figure [T71 for two redshift 
bins (0.4 < z < 0.7 and 0.6 < z < 0.9). We have not 
repeated the error estimates from the mock catalogues 
for a different mi-weighting. The {Sg\S) functions of the 
differently m^-weighted galaxy overdensity field are con- 
sistent within the la errors of the m^ = 1 weighted galaxy 
density field over the whole i5-interval probed. There is 
some indication that at positive (5, the galaxy density 
field weighted by m^ = M*^i is more biased at z 0.55, 
while the galaxy density field weighted by mi = Ls^i is 
more biased at z ~ 0.75 than the other considered types 
of galaxy density field. 

There seems to be a clear difference in the biasing re- 
sults obtained with the different mi-weights, when char- 
acterising the mean biasing function by its moments, visi- 
ble in Figure[TH| In the absence of the measured errors for 
the various mi-weightings, we have adopted the errors on 
the biasing parameters estimated for the mi = 1 weight- 
ing. They can be considered at least as an indication for 
the real errors. The b parameter increases with increasing 
z for all three types of adopted mi-weighting, hinting at a 
different speed in this evolution. The b parameters of the 
Lb- and -weighting are very similar and larger than b 
for the mi — 1 weighting at z ^ 0.75. The b parameter of 
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the i^-weighting decreases more rapidly and at z ~ 0.55 
is consistent with b parameter of the nii = 1 weighting. 
The b parameter of the M*-weighting is either larger or 
the same within the errors compared to the b parame- 
ter of the Lb- or unity-weighting. This happens due to 
higher M^/Lb ratio in high density regions progressively 
appearing as z decreases (see Bolzonclla ct al. 2009 for 
the evolution of the zCOSMOS stellar mass function in 
different environments). 

The nonlinearity parameter b/b does not evolve with 
redshift for = Lba, as we have seen already for 
rrii — 1, and it is consistent for these two galaxy den- 
sity field reconstructions. There is some evidence that 
b/b with rrii = M^^i weighting increases with decreas- 
ing 0, reaching the value of 1.02 at 2: ~ 0.55. That is 
the highest value of the nonlinearity parameter for the 
Mb < 20 — z sample and Rth = 10 h~-'^Mpc smoothing 
at all redshifts probed by the zCOSMOS. 

The behaviour of the {5g\5) function and the biasing 
parameters with the different rrii-weighting schemes for 
Mb < —20.5 — z tracer galaxies in 0.4 < z < 0.9 is 
fully consistent with the equivalent results derived with 
Mb < — 20 — z tracer galaxies (described in this section). 

6. COMPARISON TO PREVIOUS WORK 

The zCOSMOS biasing results are in qualitative agree- 
ment with the biasing framework extracted from the 
ACDM models and with a previous study of similar type 
based on the VVDS data. In order to put some tighter 
constraints on the models of galaxy formation and evolu- 
tion it is important to carry out more exact comparisons 
with previous results. For this purpose, we use only the 
linear biasing parameter, which is much easier to com- 
pare. 

6.1. Biasing parameter from the nonlinear biasing 
analysis 

In the study of the (nonline ar) biasing of the VVDS 
galaxies, iMarinoni et al.l (|2005l ) assumed a cosmological 
model described by ilm.o = 0.3, f^A^o — 0.7 and trg — 0.9. 
Moreover, they use the b parameter as the linear bias- 
ing parameter. For a proper comparison with our lin- 
ear biasing parameter 6, we infer b from the values of b 
and no nlinearity parameters published in lMarinoni et all 
(|2005l ) and correct it for the difference in the used erg nor- 
malisation, following the correction given by SBD. The 
comparison of linear biasing parameters from the zCOS- 
MOS and VVDS is presented in Figure fTUl We use the 
errors for the 6 parameter in lMarinoni et all (|2005f ) as a 
proxy for the errors of the b parameter which we have 
inferred, rescaled for the difference in erg. Even with 
this effort our comparison to the results from the VVDS 
is still only approximate. Complications arise from the 
different lumin osity complete sample s used for the analy- 
ses. Moreover. IMarinoni et al] (|2005D use a non-evolving 
magnitude to define the samples of tracer galaxies, while 
we use an evolving magnitude cut. 

As a local reference in Figure [121 we use the bias ob- 
tained fro m the bispec trum of the 2dFGRS galaxies mea- 
sured by IVerde et al. (200^. We take 61 values (fol- 
lowing the iFrv ir7}aztanagalll993l bias description, see 
Equation |4]) as a proxy for the linear biasing parameter. 



derived for L* galaxies. iNorberg et al.l ()2001[ ) detected a 
clear increase of the biasing parameter 65 (from the clus- 
tering analysis) with luminosity L of galaxies, described 

well by |l = 0.85 + 0.15-^, where 6| is the bias for the 

L* galaxies. Using this relation, we calculate the bias 
of the sample of Mb < —20 — z galaxies at the effectiv e 
redshift of 2dFGRS survey z = 0.17 (Verde et al.ll200l) . 
taking for L the median luminosity of all Mb < —20 — z 
zCOSMOS tra cer galaxies. We use the prescription from 
INorberg et al.l (2002) for the transformation from our B 
to the 2dFGRS bj filter at the median B — V colour of 
Mb < —20 — z zCOSMOS tracer galaxies with secure 
redshifts. 

It is clear that the linear bias values measured from 
the zCOSMOS and VVDS surveys at z > 0.4 are higher 
than the linear bias measured in the local universe. At 
z ~ 0.55, the mean redshift of the lowest explored bin 
of both the zCOSMOS and the VVDS surveys, there is 
an excellent agreement between the biasing parameters 
from the two surveys. However, the evolution of the lin- 
ear biasing seems to be happening at different speed at 
higher redshifts. While in zCOSMOS we detect a con- 
stant increase in the linear biasing parameter of about 
0.15 — 0.2 per 6z = 0.1 for the galaxies of the similar 
evolved luminosity, the VVDS results suggest a lower in- 
crease of the linear biasing parameter of about 0.1 or 
less per dz = 0.1. Considering the errors, the difference 
between the two results is < 2a in 0.5 < z < 1. 

The most probable explanation of this difference, apart 
from the various techniques used for the density field 
reconstruction, cosmological density parameters, galaxy 
samples etc, lies in the observed fields themselves. For 
instance, the COSMOS volume is dominated by large 
density fluctuations over the full redshift range probed 
by the zCOSMOS. Large structures domina te the whole 
zCOS MOS field even at redshifts of z ~ 0.9 (jKovac et al.l 
l2009( l. When compared to the mock catalogues, the size 
of structures in the zCOSMOS galaxy density field points 
out that the COSMOS field i s in the upper tail of the 
cosmic variance distribution (jKovac et al.ll2009D indi- 
cated also by the o ther studies (iMcCracken et al.ll2007l : 
iMeneux et aLll2OO90 . On the other hand, the VVDS field 
is in the lower tail of the cosmic variance distribution 
dM 

eneux et al.ll2009D . Moreover, as we have mentioned 
already, the field used in the VVDS analysis is smaller 
than the zCOSMOS field and the uncertainty in the red- 
shift precision in the VVDS is about 3 times lower than 
in the zCOSMOS. Nevertheless, the exact bias values and 
their redshift evolution should be explored further. 

To complete the comparison, we add in Figure fTOl linear 
bias values obtained from the clustering analysis in the 
DEEP2 survey in 0.75 < z < 1.2 with Mb < -20.77 
and Mb < —21.27 samples of galaxie s in the ra nge of 
projected scales Tp = 1 - 10 h'^Mpc (|Coil et al.i r2006). 
These bias values fall between the zCOSMOS and VVDS 
measurements, however, these are bias values obtained 
from the different statistics. We will discuss in more 
details the bias obtained from the clustering anlyses in 
the following section. 

6.2. Biasing parameter from clustering studies 

As discussed in the Introduction, clustering studies are 
commonly used to derive the linear biasing parameter. 
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One needs to keep in mind that the hnear bias inferred 
from the clustering analysis and from the biasing mo- 
ments analysis are not equivalent. In some sense, 6j is 
differential, calculated at a given distance r^, while the 

biasing parameter from our analysis b is some weighted 
average over the smoothing scale R. Moreover, in recent 
clustering studies, linear bias is derived by using simul- 
taneously correlation values at a range of r^, leading to 
a value of the linear biasing parameter more comparable 
to our approach. In addition, the difference between b 
and b^ also reflects the physical fact ors such as sto chas- 
ticity and nonlinea rity (e.g. Somcrv ille et al.ll2001h . For 
example, at z — ISomerville et al.l ( 20011 ) find that b^ 
is systematically higher than b for about 10 — 20% for 
the ACDM cosmology for the Mb < -20.27 sample of 
(mock) galaxies. 

Meneux et aT] ()2009l ) study the dependence of cluster- 
ing of the 10k zCOSMOS galaxies on their luminosity 
for the various evolving-luminosity complete samples in 
0.4 < z < 1. The measured dependence of the projected 
correlation function on the luminosity of galaxies is very 
weak and without any coherent redshift evolution in the 
amplitude or shape. From the comparison of the correla- 
tion function of Mb < —20.27— z galaxies (presented in 
iMeneux et al.l [20091 . using here h = 0.7 for magnitudes) 
to the correlation function of dark matter ([Smith et al.l 
[2003 ') in 0.5 < z < 1, these galaxies are consistent with 
biasing b^ = 1.6 at small scales, while at large scales 
R > 8 h~^Mpc the clustering analysis requires &^ ~ 2. 
The high biasing value at large scales probably reflects 
the relatively small transverse size of the zCOSMOS field 
with respect to the size of the structures, as the corre- 
la tion function is not a power law (see e.g. Figure 19 
in IMeneux eral][2009f ) . The zCOSMOS b^ inferred from 
the sample of Mb < — 20.27— z; galaxies in 0.5 < z < 1 at 



10 h-^Mpc ([Meneux et al.ll2009[ ) is larger for about 

24% and 15% than the b measured here for the samples 
of Mb < —20 — z (averaged over 0.5 < z < 0.9) and 
Mb < —20.5 — z galaxies (averaged over 0.5 < z < 1) 
with R — 10 h~^Mpc, respectively. For the smaller Vp 
scales, this b^ value is larger for about 5% and it is smaller 

for about 6% than b for the same zCOSMOS samples as 
above. 

It is of interest to note here that IMeneux et al. (2009) 
compare the cor relation function of t he zCOSMOS and 
VVDS galaxies ([Meneux et al.ll2008f) . finding that the 
bias inferred from clustering of log(M,/M0) > 10 galax- 
ies in 0.5 < z < 1 is systematically higher for the zCOS- 
MOS than VVDS galaxies. The observed difference is 
fully consistent with the difference between the linear 
biasing parameters derived from the nonlinear biasing 
analysis in these two surveys, as shown here (e.g. Fig- 
ure [T9|). 

Using the full DEEP2 sample, ICoil et al.l ([2006") mea- 
sure the increase of 6^ with luminosity on both small and 
large scales, where the trend is stronger on smaller scales. 
The measured bias b^ takes values from 1.42 ± 0.04 for 
the Mb < -19.77 sample at 0.75 < z < 1 to 1.67 ± 0.06 
for the Mb < -20.77 sample at 0.75 < z < 1.2 or to 
1.74 ± 0.05 for the Mb < -21.27 sample of galaxies at 
0.75 < z < 1.2 (using here h = 0.7 for magnitudes), de- 
rived simultaneously for the range of scales = 1 — 10 



h -'^Mpc. We h ave increased the bias values and their 
errors quoted in lCoil et all ([2006f ) by - 13% in order to 
correct for the different erg u sed. The va l ues o f the most 
luminous DEEP2 samples of ICoil et all ([2006f) are com- 
parable to our b values of the two most luminous samples 
at z > 0.6, and they are higher than any of the b values 
inferred from the VVDS sample of Mb < —20.77 galaxies 
in 0.7 < z < 1.5 (see Figure [Tg]). Moreover, the DEEP2 
sample is complete for the red A/r < —21.27 galaxies 
only up to z = 1.05 ([Coil et al.ll2008D . T hese galaxies are 
more biased than the blue galaxies (e.g. lCoil et al.|[2008[) 
and therefore the value of bias for the Mb < —21.27 
DEEP2 galaxies can be partially lower due to this effect 
(our luminosity complete samples are chosen to be also 
colour complete). 

The comparison of the zCOSMOS linear biasing mea- 
sured in this work to the previous biasing analyses based 
both on the clustering and moments statistics, leads us 
to conclude that the cosmic variance is the main contrib- 
utor to the different bias values and the different rate 
of the redshift evolution of bias of galaxies with similar 
luminosities in the existing 0.4 ^ -z ^ 1-5 spectroscopic 
surveys. 

6.3. Bias of the dark matter haloes 

From the observed bias of the zCOSMOS galaxies we 
can infer a characteristic mass of dark matter haloes that 
host these galaxies. Under the assumption that the ob- 
served galaxies are central galaxies of a halo, the mea- 
sured bias (which as we showed does not depend on 
the scale of 8 to 12 h~^Mpc used to reconstruct the 
galaxy density field) can be matched to the bias of the 
dark matter haloes. With this assumption (no satellite 
galaxies), the infered characteristic halo mass at a given 
bias value will be higher than the true halo mass (e.g. 
[Zheng et al.l2007( ). In the zCOSMOS survey, the fraction 
otiAB < 22.5 galaxies in groups is increasing from ~ 15% 
at z - 0.8 to 25% at z - 0.1 ([Knobel et al.ll2009f ). and 
this fraction is strongly dependent on galaxy properties 
(e.g. stellar mass, see Figure 12 in Kovac et al. 2009). 
The additional uncertainty comes from the fact that we 
are comparing bias of galaxies above a given luminosity 
with a bias of haloes of a given mass. Having in mind all 
the approximations mentioned above, the results on the 
typical halo masses of the zCOSMOS galaxies should be 
understood only as indicative. 

We use two approximations for the bias of haloes of a 
give n mass Mh at r e dshift z given bv lSheth et aH ([200I 
and IPillepich et al.l ([2008D , to also highlight theoretical 
uncertainties. First, we calculate for the adopted cosmol- 
ogy the halo bias for the range of halo masses at redshifts 
0.55, 0.65, 0.75 and 0.85, which are the mean redshifts of 
the redshift intervals used in our analysis. We start from 
the simplest scenario, in which galaxies of a given type 
reside always in haloes of the same mass. For the sample 
of Mb < —20 — z galaxies and R — 10 h^^Mpc, we find 
by minimisation of differences between the bias of galax- 
ies and haloes at different redshifts that the best fit mass 
of dark haloes to host this sample is ^ 3 x 10^^ h~^Mp, 
or ^ 6 X 10^^ h-^Mfn, w hen using the lSheth efall ([200 11 1 
or iPillepich et al] ((2D08) bias expressions, respectively. 

The redshift evolution of linear biasing of ~ 3 x 
10^2 h-^Mg and - 6 X lO^^ h-^M© haloes following 
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iSheth et all (pOOl and iPillepich et all (iool . respec- 
tively, is presented in Figure[20l along with the measured 
redshift evolution of linear bias of the Mb < —20 — z 
zCOSMOS galaxies. The models of halo biasing approx- 
imately describe the observed evolution of galaxian bias- 
ing, where the difference falls almost completely within 
the 1(7 cosmic variance errors of the galaxy bias. How- 
ever, our result indicates that the evolution of galaxy bias 
is faster than the evolution of halo bias for the haloes of 
the consi dered mass as p redicted by the both models, 
with Pil lepich et al.l ()2008[ ) bias evolution matching bet- 
ter the observed trend. 

The considered (halo) biasing models include the ef- 
fect of the merging of haloes. For a comparison, we cal- 
culate also the evolution of biasing vising the so called 
"galaxy conserving" model (e.g. Fry 1996), in which the 
number of galaxies is preserved over cosmic time (no 
merging). Here, we assume that the model bias at the 
redshifts of zCOSMOS observations is given by the bias 
of Mb < -20- z zCOSMOS galaxies at the i? = 10 
h~^Mpc scale. The results are included in Figure 
The conserving model produces too high values of bi- 
asing at low z. This is a known result in the biasing 
analysis, indicating that merging is an important ingre- 
dient in the biasing of cosmic structures. Commonly, 
the difference in the evolution of halo and galaxy biasing 
(e.g. black points and red/blue curves in Figure [^0)1 are 
attributed to the different timescales of mergers of galax- 
ies and haloes, as well as the evolution in mass-to-light 
ratios between haloes and galaxies (e.g. Somcrville et al.l 
[2001). For the precise answer on the difference in the 
evolution of galaxian and halo bias, the broader baseline 
in redshift is needed. 

If we consider a different scenario, in which the galax- 
ies defined by their evolving _B-band luminosity reside in 
haloes of different mass at different redshifts, we find that 
the characteristic mass of a halo (i.e. the halo mass at 
which the halo bias matches the bias of galaxies) to host a 
Mb < — 20 — z galaxy increases from ~ 1.8 x 10^^ h~^M0 
at z ~ 0.55 to 4.3 x 10^^ h~^ M(7) at z ~ 0.75 whe n using 
the halo bias approximation of lSheth et al.l (|2001l contin- 
uous curves in Figure [2T|) or from ~ 4.3 x 10^^ h~^MQ at 
z - 0.55 to 7.7 X 10^2 h~^ M(7i at z ~ 0.75 when u sing the 
halo bias approximation of IPillepich et al.l ()2008l . dashed 
curves in Figure I^T|) . This result indicates that galaxies 
defined by the same — Az evolved luminosity reside on 
average in more massive haloes at higher z (where z < 1). 
We do not include errors in the quoted characteristic halo 
masses, and therefore these values should not be used to 
infer the exact rate in evolution in the halo mass-to-light 
ratio, particularly having in mind the discrepancy from 
the zero mean overdensity field in the 0.4 < z < 0.7 bin, 
which can largely affect the implied evolution. Physi- 
cally, the inferred trend in the halo mass-to-light ratio is 
consistent with a scenario in which the star formation, 
as traced by the B-band luminosity, shifts from more 
massive to less massive haloes with decreasing redshift. 
This is similar to the evolution in halo mas s-to-light ra- 
tio fo und between t he DEEP2 and SPS S (|Zheng et al.1 
I2007f ). Specifically, IZheng et al.1 (|2007D find that the 
mean luminosity of the central galaxy increases with halo 
mass at both redshifts, and the central galaxies reside 
in the haloes a few times more massive at z '--^ 1 than at 
z 0. However, the physical interpretation of the results 



above is hampered by the fact that we fand lZheng et al.1 
l2007f ) have used the luminosity complete samples, defined 
here by their evolving _B-band luminosity, which is very 
sensitive to the recent star formation history. Preferen- 
tially, one should use stellar masses to define complete 
samples. Moreover, as seen in Section 15.51 the biasing 
inferred from the stellar mass-weighted density field is 
higher, and its evolution can be well described by a bias 
of the halo of a constant characteristic mass. The con- 
nection between galaxies and haloes inferred from the 
differently weighted density field has not been explored 
yet. 

From Figure it is also noticeable that the models for 
the halo bias predict higher biasing of haloes of higher 
mass. Moreover, the difference between the biasing of 
haloes of different masses increases with redshift, reflect- 
ing the fact that at high redshift more massive haloes are 
formed in the higher, and more rare density peaks, and 
therefore they are more biased tracers of the underlying 
mass distribution at higher redshifts. We measure higher 
linear biasing for more luminous galaxies, and therefore 
the characteristic halo mass of more luminous galaxies is 
higher. For examp le, when using the bias expression by 
ISheth et al.l (|2001h . at z ~ 0.55 the characteristic mass 
of haloes which host Mb < 19.5 — z and Mb < —20.5 — z 
zCOSMOS galaxies is 1.3 x lO^^ Yi'^Mq and 3.5 x lO^^ 
h-^Mg, respectively. At z ~ 0.75, Mb < -20.5 - z 
zCOSMOS galaxies are hosted by dark matter haloes of 
mass o f 6.9 x 10^^ h~^M(7). W hen using the bias expres- 
sion bv IPillepich et all (120081 ). the corresponding charac- 
terstic masses are 3.5 x 10^^ h'^Mg, 7.1 x 10^^ h-^Mg 
and 1.1 X 10^3 i^-^Mq for Mb < 19.5 - z at z - 0.55, 
Mb < -20.5 - z at z ~ 0.55 and Mb < -20.5 - z at 
z ^ 0.75 zCOSMOS galaxies, respectively. 

7. SUMMARY AND CONCLUSIONS 

In this work, we make use of the reconst ructed over- 
densi ty field in the zCOSMOS volume (see iK ovac et al.l 
l2009f ) to derive the conditional mean function {Sg\S) = 
b{S, z, R)6 and the second moments of the mean bias- 
ing function b{S,z,R). For this purpose we employ the 
density field on a grid reconstructed by using the three 
dimensional distances between galaxies and grid points 
and counting the objects within a spherical to p-hat aper- 
ture. We implement a novel method ZADE ((Kovac et al.l 
|2009| ) to account for galaxies not yet observed spectro- 
scopically in the selected samples of galaxies used to re- 
construct the density field. For a biasing analysis, the 
main advantage of the ZADE method is that in a statis- 
tical sense, the mean intergalaxy separation is that of all 
galaxies in the selected galaxy sample, and not only of 
the sample of galaxies with spectroscopic redshifts. 

We have carried out a number of tests on the mock cat- 
alogues to assess various errors which are going to affect 
our biasing analysis. Particularly, we have empirically es- 
timated uncertainties due to cosmic variance, shot noise 
errors and the density field reconstruction errors. 

• Cosmic variance errors cause a spread in the {Sg\S) 
function: for a given S there is a range of dg values 
measured in the mock catalogues. Quantifying this 
spread by the standard deviation cr of log(l -I- Sg), 
we find that a is largest in the most underdense 
regions where a ^ 0.1, becomes lower at the inter- 
mediate S values, a < 0.05, and increaseas again 
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in the most overdense regions up to cr ~ 0.05 at a 
given S. 

• The shot noise (discrete galaxy samphng) errors 
modify significantly the shape of the {Sg\5) function 
in the most underdense regions, making the local 
bias b{S, z, R) values in the same regions to appear 
higher. 

• Reconstruction errors are relevant only in the un- 
derdense regions. At the reconstructed value of 
log(l + 5g) = — 1 in the galaxy density field, the re- 
construction error can cause an uncertainty of the 
order of 0.1 in the matter density field log(l + S). 

• The b parameter increases due to the shot noise and 
reconstruction errors. The b/b parameter is not 
susceptible to either of these errors. The cosmic 
variance causes a spread in the measured values of 
both of these parameters. 

We can summarise our main findings in the biasing 
analysis of the 10k zCOSMOS galaxies as follows: 

• The conditional mean function {Sg\6) has a charac- 
teristic shape as described below. In most under- 
dense regions, the mean biasing function vanishes. 
At some 6 < 0, the mean biasing function appears 
and then rises sharply in the underdense regions, 
with the local slope of the biasing function larger 
than unity. Starting from around mean density 
and towards higher overdensities, the {6g\5) func- 
tion closely follows a linear relation Sg = b6 with b a 
constant. In the most overdense regions zCOSMOS 
galaxies are antibiased, i.e. locally b{d,z,R) < 1. 
This is true for all samples of tracer galaxies used. 
The conditional mean function {Sg\S) is clearly non- 
linear in the most overdense and underdense re- 
gions. 

• There is a detectable change in the shape of the 
{6g\S) function with redshift in the overdense re- 
gions. For a given population, galaxies become 
more biased tracers of the matter in the regions 
of mean and mildly positive overdensities as red- 
shift increases from 0.4 to 1. There is an indication 
of an evolution in the value of S in the overdense 
regions, at which galaxies become antibiased. For 
a given population of tracer galaxies, this happens 
at higher S for higher z. Taking into account all 
the sources of errors, we cannot discriminate if the 
shape of the {Sg\S) function in the underdense re- 
gions stays constant or undergoes some evolution 
in 0.4 < 2; < 1. It is worth stressing that redshift 
evolution may be not visible due to the overlapping 
redshift bins in the analysis. 

• The {dg\d) function shows some dependence on 
the scale in the most overdense regions, more pro- 
nounced at high redshift, such that at smaller scales 
galaxies are more biased tracers of the underlying 
matter distribution. However, this is not highly 
statistically significant. 



• When comparing the {Sg\S) function for the differ- 
ent luminosity selected tracer galaxies, there is an 
indication that more luminous galaxies are more 
biased than less luminous galaxies in the regions 
from about mean density to the highly overdense 
regions. 

• The linear bias between the 10k zCOSMOS galax- 
ies and mass is increasing with redshift. There 
is some evidence that more luminous galaxies are 
more biased tracers of matter. We do not detect 
any significant dependence of the linear biasing pa- 
rameter on the scale at which we measure galaxy 
overdensity fields. The nonlinearity parameter is 
of the order of a few percent. It does not change 
with the redshift, with the smoothing scale or with 
the luminosity of the tracer galaxies. 

• The linear bias of the stellar mass-weighted den- 
sity field is either larger or the same within the 
errors compared to the linear bias of the -B-band 
luminosity- or unity-weighted density field. 

• By comparing the galaxy biasing to the halo bi- 
asing, using the appro ximation for h a lo bia s of 
ISheth et al.l (|2001[ ) and Pi Uepich et all (|2008f) . we 
infer that the Mb < -20 - z zCOSMOS galax- 
ies in 0.4 < z < 1 reside in dark matter haloes 
with a characteristic mass of '--^ 3 or ~ 6 x 10^^ 
h~^MQ, using these two models, respectively. One 
would need to work with the stellar mass complete 
samples to obtain the evolution in characteristic 
halo mass whose physical interpretation is not am- 
biguos. 

Broadly speaking, our results are in line with find- 
ings from the previous stud y of the nonlinea r biasing of 
high redshift galaxies (Mar inoni et al.ll2005l I2OO8;) and, 
qualitatively, they follow th e biasing history outlined 
by the theoretical wor ks (e.g. iBlanton et al.lfl999l l2000l : 
ISomerville et al.lf2001f ). When going into details, there 
are a number of discrepancies which neeed to be solved, 
such as the exact "speed" in the evolution of the linear 
biasing parameter or the dependence of biasing on the 
luminosity of tracer galaxies. While the current biasing 
results from the 2: < 1.5 surveys are important as they 
support the general picture of the biased galaxy forma- 
tion and provide a framework for the future theoretical 
work, the fine tuning of the galaxy formation models is 
still hampered by the limitations of the existing spec- 
troscopic surveys at these redshifts, particularly by the 
cosmic variance. 
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TABLE 1 

Summary of the measured biasing parameters. 




Fig. 1. — Unity-weighted (m^ = 1) mean volume density of the zCOSMOS tracer galaxies in 0.1 < z < 1. The dotted lines are the 
mean volume densities of the zCOSMOS tracer galaxies obtained by dividing the number of galaxies with volume corresponding to redshift 
bins of 0.05. The continuous lines are smoothed mean densities, obtained by adding AV{z)/Vmax contributions of each tracer galaxy in 
redshift bins 0.002 wide. The magenta, blue, red and green curves are for the flux limited, and Mb < —19.5 — 2, Mb < —20 — z and 
AIb < —20.5 — 2 luminosity complete samples, respectively. See text for more details. 
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Fig. 2. — Example of the PDF of the galaxy density contrast field using a mock catalogue. The density field is reconstructed using the top- 
hat smoothing filter of 5, 8 and fO h~^Mpc going from the top to the bottom panel. The black histogram corresponds to the reconstruction 
with the 40k catalogue (every galaxy has a measured spectroscopic redshift) and the red histogram corresponds to the reconstruction with 
the 10k+30kZADE catalogue (lOk-like sample of galaxies with a measured spectroscopic redshift, the rest of the Ij^g < 22.5 galaxies have 
the ZADE-modified photometric redshift). Binning is carried out in log{l + Sg) units. At the current status of the zCOSMOS survey, we 
need scales of at least 8 h~^Mpc to reconstruct the density field for the biasing analysis with acceptable errors at every Sg up to 2 ~ 1. 
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Fig. 3. — Effect of the cosmic variance errors on the conditional mean function {5g\5). At every S value, cosmic variance is causing a 
spread in the corresponding Sg values. The mean {5g\5) functions are obtained from the 12 mock catalogues, plotted as the thin dotted 
lines, where the mock galaxy density field has been reconstructed with the fiux limited sample of galaxies and top-hat filter of 10 h~^Mpc 
in 0.4 < z < 0.7. In all figures containing the {5g\5) function, starting with this, we mark the case of a no-biasing bi^ = 1 with the black 
dotted line. The cross in the middle of the panels is drawn for a reference and marks the Sg = S = case. The standard deviation of 5g 
values in the mocks at every S are plotted with the thick continuous lines, centred at the mean 5g. This is the effective cosmic variance 
noise expected in a single reconstruction, i.e. in the actual data. 
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Fig. 4. — Effect of galaxy sampling on the conditional mean function (<5g|<5). The (<5g|(5) function between galaxies (Ms < —18) and 
mass for a top-hat smoothing of 8 h^^Mpc in two redshift bins: Az = 0.4 — 0.7 and Az = 0.7 — 1.0 is presented in the left and right 
panel, respectively. The (5g|<5) function with the full sampling in the given mock catalogue {I ^ 2.7 h~^Mpc) is presented with the black 
continuous line. The various curves correspond to the different mean galaxy separation I of 3 (yellow), 4 (pink), 5 (red), 6 (violet), 7 (blue) 
and 8 (cyan) h~^Mpc. 
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Fig. 5. — Summary of the effect of the sampling of galaxies on the second moments of the mean biasing function. Individual points 
correspond to the mean biasing parameters; linear bias h (top panels) and nonlinearity (bottom panels) obtained by averaging results 
from 12 mocks, plotted as a function of the mean intergalaxy separation I. The vertical bars correspond to the standard deviation of the 
parameters from these 12 mocks. The colour coding is the same as in the previous figure. The four panel plots on the left and right side 
are obtained for the galaxy density field reconstructed on i? = 8 h~^Mpc and i? = 10 h~^Mpc, respectively. In each of the four panels, 
the left hand plots refer to 0.4 < z < 0.7 and the right hand plots refer to 0.7 < z < 1. 
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Fig. 6. — Reconstruction errors in the mean biasing analysis for the fiux complete samples. Lower panels: the {Sg\5) function of the 
overdensity field of galaxies {Iab < 22.5) for a top-hat smoothing of 8 h~^Mpc for the lOk+SOkZADE reconstruction (violet) and 40k 
reconstruction (black). The curves are obtained by averaging results from 12 mock catalogues of the same type. Upper panels: the scatter 
in the corresponding (i5g|<5) functions plotted below. The scatter at fixed 5 values is calculated as the standard deviation (cr) of 5g values 
from 12 40k- type mocks (black) and from 12 10k+30kZADE-type mocks (violet). 
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Fig. 7. — Reconstruction errors in the mean biasing analysis for the luminosity complete samples. The curves have the same meaning as 
in Figure |S] but the zCOSMOS and mock overdensity field is reconstructed with the Mg < —20 — z galaxies. 



- •40k 

1.4 - • 10k+30kZADE 



1.2 



1 - 
1.04 r 



I I I I I I 



I,B<22.5 
R,„=8h-'Mpo 



I I I I I 



,1.02 r 



il il II .1 



0.4 0,5 



0.8 0,9 



- •40k 

1,4 - • 10k+30kZADE 



1,2 



1 - 
1,04 r 

,1,02 r 



I I I I I I 



Mj<-20-z 
R,„=8h-'Mpc 



I I I I I 



0,4 0,5 0,6 0,7 



0,8 0,9 



Fig. 8. — Summary of the effect of the reconstruction errors on the second moments of the mean biasing function. Individual points 
correspond to the mean biasing parameters: b (top panels) and b/b (bottom panels) obtained by averaging results from 12 mocks. The 
vertical bars are the corresponding standard deviations. 
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Fig. 9. — zCOSMOS overdensity field reconstructed with R = S h~^Mpc and flux limited tracer galaxies in 0.4 < z < I. The colour 
scale on the bottom is given in the 1 + S units. The horizontal axis are RA and DEC, the vertical axis is redshift. The size of the box is 
35, 40, 45 and 50 h~^Mpc along RA and DEC axis from lower to higher redshift and ~ 0.15 along redshift. We plot only structures above 
the mean density (1 + 5>1), in order to increase the visibility. 
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Fig. 10. — Conditional mean function {Sg\S) estimated for different eg values. The parameter erg increases from 0.7 to 1 in steps of 0.05. 
The resulting {Sg\5) function is presented in blue for erg = 0.7, cyan for eg = 0.75, black for erg = 0.8, green for erg = 0.85, orange for 
erg = 0.9, red for erg = 0.95 and violet for erg = 1. 
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Fig. 11. — Conditional mean function {Sg\5) for the density field of the zCOSMOS galaxies {Mb < —20 — z) obtained by smoothing on 
scales of 8 h~^Mpc. The red curve corresponds to the linear biasing case Sg = bS. The different panels are for the different redshift bins: 
0.4 < 2 < 0.7, 0.5 < z < 0.8 and 0.6 < 2 < 0.9 from the left to the right, respectively. 
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Fig. 13. — Conditional mean function {Sg\5) for the density field of the zCOSMOS galaxies (Mb < —20 — z) obtained by smoothing on 
various scales. The resulting functions are presented with the black continuous, cyan dotted and violet dashed lines for the scale of 8, 10 
and 12 h~^Mpc, respectively, in two redshift intervals; 0.4 < z < 0.7 (left) and 0.6 < z < 0.9. 
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Fig. 14. — Biasing parameters for the 10k zCOSMOS overdensity field {Mg < —20 — z) calculated for the various smoothing scales. The 
redshift evolution of the linear biasing parameter b is presented in the top panel and the redshift evolution of the nonlinearity parameter 
b/b is presented in the bottom panel. The black, cyan and violet symbols represent the parameters for the scale of 8, 10 and 12 h~^Mpc, 
respectively. The errors are calculated as the standard deviation of the corresponding biasing parameters from 12 40k-type mocks. Some 
points and their errors are displaced along redshift-axis from the mean redshift in the bin where analysis was carried out for the sake of 
clarity. 
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Fig. 15. — Conditional mean function (Sg\5) for the density field of the zCOSMOS galaxies of various luminosity thresholds obtained by 
smoothing on a scale of 10 h~^Mpc. The resulting functions are presented with the black continuous, yellow dotted, red dashed and green 
dot-dashed lines for the samples of Ij^g < 22.5, Mg < —19.5 — z, Mg < —20 — z and Mg — 20.5 — z galaxies, respectively, in two redshift 
intervals: 0.4 < z < 0.7 (left) and 0.6 < z < 0.9 (right). For the luminosity complete samples, only results for the samples for which we are 
complete in a given redshift interval are presented. 
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Fig. 16. — Biasing parameters for the 10k zCOSMOS overdensity field of the various luminosity thresholds. The black, yellow, red and 
green symbols represent the parameters for the samples of I^g < 22.5, Mg < —19.5 — z, Mg < —20 — z and Mg — 20.5 — z galaxies, 
respectively. For the luminosity complete samples, only results for the samples for which wo arc complete in a given redshift interval are 
presented. Details are as in Figure [141 
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Fig. 17. — Conditional mean function {5g\5) for the density field of the zCOSMOS galaxies {Mb < —20 — 2) obtained by smoothing on a 
scale of 10 h~^Mpc with the different weighting schemes nii (see Equation [Sj for the density field reconstruction. The resulting functions 
are presented with the black continuous, green dotted and magenta dashed lines for the mi = 1, rm = Lg i and rm = M*^^ weighting 
schemes, respectively, in two redshift intervals: 0.4 < z < 0.7 (left) and 0.6 < z < 0.9 (right). 
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Fig. 18. — Biasing parameters for the 10k zCOSMOS overdensity field {Mb < -20- z and i? = 10 h~lMpc) calculated with the different 
weighting schemes nii (see Equation [Sj for the density field reconstruction. The black circles, green traingles an d p ink squares represent 
the parameters for the samples of m,i = 1, nii = Lg^i and rrii = M^^i weighting schemes. Details are as in Figure [HI 
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Fig. 19. — Comparison of the zCOSMOS linear biasing parameters to the bias values from similar analyses available in the literature. 
The filled squares correspond to the zCOSMOS measurements of b, the empty circles are b values based on the nonlinear biasing analysis in 
the VVDS (Marinoni et al. 2005) and the cross is the bi value from the 2dFGRS (Verde et al. 2002). The zCOSMOS points are presented 
at the mean redshift of the bins Az = 0.3 wide. The three lower z points are calculated for the Mg < —20 — z sample, while the z ~ 0.85 
point is for the Mb < —20.5 — z sample and the smoothing scale is = 10 h~^Mpc for all the points. The biasing values from the VVDS 
are inferred from the overdensity field reconstructed using a sample with Mg < —20.77 on i? = 5 h~^Mpc scale in 0.4 < z < 0.7, and on 
R = 10 h-^Mpc scale in 0.7 < 2; < 0.9, 0.9 < 2 < 1.1, 1.1 < z < 1.3 and 1.3 < 2 < 1.5. The VVDS points are plotted at the centre of 
the corresponding bin, with the exception of the lowest z point, offset along the redshift axis due to clarity by —0.05. The 61 value from 
the 2dFGRS is recalculated for Mg < —20 — 2 galaxies at 2 = 0.17, whic h is the effective redshift of 2dFGRS. For a comparison, we add 
linear bias values obtained from the clustering statistics in the DEEP2 llCoil et al.ll200(l) . represented as triangles. The DEEP2 bias is 
plotted at the mean redshift of the 0.75 < 2 < 1.2 interval used for the analysis, offset along the redshift axis due to clarity by —0.02 for 
the Mg < —20.77 sample and by +0.02 for the Mg < —21.27 sample of galaxies. See text for more details. 




Fig. 20. — Redshift evolution of the bias of Mg < —20 — 2 galaxies (filled circles), 3 X lO^'^ h ^Mq dark matter haloes following the 
IS heth et al] (2001. S2001) expression for halo biasing (continuous line) and 6 X 10^^ h ^ Mq dark matter haloes following the Piilepi ch et al.l 

f008, P2009) expression for halo biasing (dashed line). The evolution of biasing of haloes using the so called "galaxy conserving" model 
rv 1996, F1996) is calculated assuming that at the redshift of zCOSMOS observations the model bias has the same value as the bias of 
Mg < -20 - z zCOSMOS galaxies at the _R = 10 h'^Mpc scale (dotted fines). 
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Fig. 21. — Bias of the dark matter haloes at reds hifts used in t his an alysis. The c ontinuous and dashe d curves correspond to the bias of 
dark matter haloes of a given halo mass as given bv lSheth et~all ll200ll . S2001) and Pillcpi ch et al.l l l200g . P2009), respectively. The results 
are calculated for redshifts z = 0.55, z = 0.65, z = 0.75 and z = 0.85 (mean redshifts of the bins used in the presented biasing analysis) for 
the curves presented from the bottom to the top, as indicated by the arrow along the given redshift values. The filled symbols (squares and 
triangles) are b values from the zCOSMOS analysis with the Mb < —20 — z sample on ij = 10 h~^M pc scale, plotted at the values of a 
dark matter halo mass of the same bias as galaxies at the corresponding redshift. The squares are for the lSheth et al.l (120011 ) approximation 
and the triangles are for the , Pillcpich et al.. (.2008.1 approximation. 



